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O ' Abstract 
(N 

There has been much activity in the measurement of the elastic electromagnetic proton and neutron 
form factors in the last decade, and the quality of the data has been greatly improved by performing double 
polarization experiments, in comparison with previous unpolarized data. Here we review the experimental 
data base in view of the new results for the proton, and neutron, obtained at MIT-Bates, MAMI, and JLab. 



The rapid evolution of phenomenological models triggered by these high-precision experiments will be 
discussed, including the recent progress in the determination of the valence quark generalized partem distri- 
butions of the nucleon, as well as the steady rate of improvements made in the lattice QCD calculations. 
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1 Introduction 



The characterization of the structure of the nucleon is a defining problem of hadronic physics, much like the 
hydrogen atom is to atomic physics. Elastic nucleon form factors (FFs) are key ingredients of this character- 
ization. As such, a full and detailed quantitative understanding of the internal structure of the nucleon is a 
necessary precursor to extending our understanding of hadronic physics. 

The electromagnetic (e.m.) interaction provides a unique tool to investigate the internal structure of the 
nucleon. The measurements of e.m. FFs in elastic as well as inelastic scattering, and the measurements of 
structure functions in deep inelastic scattering of electrons, have been a rich source of information on the 
structure of the nucleon. 

The investigation of the spatial distributions of the charge and magnetism carried by nuclei started in the 
early nineteen fifties; it was profoundly affected by the original work of one of its earliest pioneers, Hofstadter 
and his team of researchers [Hof53bJ, at the Stanford University High Energy Physics Laboratory. Quite early 
the interest turned to the nucleon; the first FF measurements of the proton were reported in 1955 [Hof55], and 
the first measurement of the neutron magnetic FF was reported by Yearian and Hofstadter IYea581l in 1958. 
Simultaneously much theoretical work was expanded to the development of models of the nucleus, as well as 
the interaction of the electromagnetic probe with nuclei and the nucleon. The prevailing model of the proton 
at the time, was developed by Rosenbluth HRos50H , and consisted of a neutral baryonic core, surrounded by a 
positively charged pion cloud. 

Following the early results obtained at the Stanford University High Energy Physics Laboratory, simi- 
lar programs started at several new facilities, including the Laboratoire de l'Accelerateur Lineaire in Orsay, 
(France), the Cambridge Electron Accelerator, the Electron-Synchrotron at Bonn, the Stanford Linear Acceler- 
ator Center (SLAC), Deutsches Elektronen-Synchrotron (DESY) in Hamburg, the 300 MeV linear accelerator 
at Mainz, the electron accelerators at CEA-Saclay, and at Nationaal Instituut voor Kernfysica en Hoge Energie 
Fysica (NIKHEF). The number of electron accelerators and laboratories, and the beam quality, grew steadily, 
reflecting the increasing interest of the physics problems investigated and results obtained using electron scat- 
tering. The most recent generation of electron accelerators, which combine high current with high polarization 
electron beams, at MIT-Bates, the Mainz Microtron (MAMI), and the Continuous Electron Beam Accelerator 
Facility (CEBAF) of the Jefferson Lab (JLab), have made it possible to investigate the internal structure of 
the nucleon with unprecedented precision. The CEBAF accelerator adds the unique feature of high energy 
which allows to perform measurements of nucleon e.m. FFs to large momentum transfers. Sizable parts of 
the programs at these facilities were and are oriented around efforts to characterize the spatial distribution of 
charge and magnetization in nuclei and in the nucleon. 

The recent and unexpected results from JLab of using the polarization transfer technique to measure the 
proton electric over magnetic FF ratio, Ge p /Gm p HJonOOl |Gay02[ |Pun05L has been the revelation that the 
FFs obtained using the polarization and Rosenbluth cross section separation methods, were incompatible with 
each other, starting around Q 2 = 3 GeV 2 . The FFs obtained from cross section data had suggested that 
Ge p ~ Gm p I 'fJ> P , where [L v is the proton magnetic moment; the results obtained from recoil polarization data 
clearly show that the ratio Ge p /Gm p decreases linearly with increasing momentum transfer Q 2 . The numerous 
attempts to explain the difference in terms of radiative corrections which affect the results of the Rosenbluth 
separation method very significantly, but polarization results only minimally, have led to the previously ne- 
glected calculation of two hard photon exchange with both photons sharing the momentum transfer. 

These striking results for the proton e.m. FF ratio as well as high precision measurements of the neutron 
electric FF, obtained through double polarization experiments, have put the field of nucleon elastic e.m. FFs 
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into the limelight, giving it a new life. Since the publication of the JLab ratio measurements, there have 
been two review papers on the subject of nucleon e.m. FFs [Gao03, Hyd04|, with a third one just recently 
completed [Arr06J. The present review complements the previous ones by bringing the experimental situation 
up-to-date, and gives an overview of the latest theoretical developments to understand the nucleon e.m. FFs 
from the underlying theory of the strong interactions, Quantum Chromodynamics (QCD). We will focus in this 
review on the space-like nucleon e.m. FFs, as they have been studied in much more detail both experimentally 
and theoretically than their time-like counterparts [Bal05|. We will also not discuss the strangeness FFs of 
the nucleon which have been addressed in recent years through dedicated parity violating electron scattering 
experiments. For a recent review of the field of parity violating electron scattering and strangeness FFs, see 
e.g. Ref. HBei05L 

This review is organized as follows. Section 2 is dedicated to a description of the beginning of the field of 
electron scattering on the nucleon, and the development of the theoretical tools and understanding required to 
obtain the fundamental FFs. Elastic differential cross section data lend themselves to the separation of the two 
e.m. FFs of proton and neutron by the Rosenbluth, or LT-separation method. All experimental results obtained 
in this way are shown and discussed. 

Section 3 discusses the development of another method, based on double polarization, either measuring 
the proton recoil polarization in ep — > ep, or the asymmetry in ep — > ep. The now well documented and 
abundantly discussed difference in the FF results obtained by Rosenbluth separation on the one hand, and 
double polarization experiments on the other hand, is examined in section 13.41 The radiative corrections, 
including two-photon exchange corrections, essential to obtain the Born approximation FFs, are discussed in 
details in section [331 

In Section 4, we present an overview of the theoretical understanding of the nucleon e.m. FFs. In Sect. 4.1, 
we firstly discuss vector meson dominance models and the latest dispersion relation fits. To arrive at an 
understanding of the nucleon e.m. FFs in terms of quark and gluon degrees of freedom, we next examine in 
Sect. 4.2 constituent quark models. We discuss the role of relativity when trying to arrive at a microscopic 
description of nucleon FFs based on quark degrees of freedom in the few GeV 2 region. The present limitations 
in such models will also be addressed. In Sect. 4.3, we highlight the spatial information which can be extracted 
from the nucleon e.m. FFs, the role of the pion cloud, and the issue of shape of hadrons. Sect. 4.4 discusses the 
chiral effective field theory of QCD and their predictions for the nucleon e.m. FFs at low momentum transfers. 
Sect. 4.5 examines the ab initio calculations of nucleon e.m. FFs using lattice QCD. We will compare the most 
recent results and the open issues in this field. We also explain how the chiral effective field theory can be 
useful in extrapolating lattice QCD calculations for FFs, performed at larger than physical pion mass values, to 
the physical pion mass. In Sect. 4.6, we present the quark structure of the nucleon and discuss how the nucleon 
e.m. FFs are obtained through sum rules from underlying generalized (valence) quark distributions. We show 
the present information on GPDs, as obtained from fits of their first moments to the recent precise FF data set. 
Finally, in Sect. 4.7, we outline the predictions made by perturbative QCD at very large momentum transfers 
and confront them with the FF data at the largest available Q 2 values. 

We end this review in Section 5 with our conclusions and spell out some open issues and challenges in this 
field. 
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2 Nucleon form factors from eN cross sections 



In this section we outline the development of what was, in the early nineteen fifties, a new and exciting field of 
investigation of the structure of nuclei, using the elastic scattering of electrons with several hundreds of MeV 
energy. We also discuss the evolution of the Rosenbluth separation method to its present form, and show all 
FF results obtained using this method for both the proton and the neutron. 

In the late forties, several papers had pointed out the possibility of measuring the shape and size of nuclei 
by observing deviations from Mott scattering by a point charge; most influential were the papers by Rose 
[Rose47], who argued that "high energy" electrons would be most suited for such studies, with 50 MeV a best 
value; and by Rosenbluth HRos50l for the proton, who provided explicit scattering formula taking into account 
both charge and the anomalous magnetic moment, with the use of "effective" charge and magnetic moment. 

An early report of work done at the Stanford University High Energy Physics Laboratory at energies larger 
than 100 MeV, was reported by Hofstadter, Fechter and Mclntyre [Hof53a], who detected deviations from 
scattering by a point charge in carbon and gold. The first review paper of the field, written by Hofstadter in 
1956 [Hof56J included measurement of the proton FF, up to a momentum transfer squared of q 2 = 13.3 fm~ 2 , 
or 0.52 GeV 2 . 



2.1 Early nucleon structure investigations 



In the middle nineteen fifties, it had been known for more than 20 years that the proton could not be just a 
mathematical point charge and point magnetic moment. Indeed the measurement of the proton's magnetic 
moment by Stern HSte33l had revealed a value ~2.8 times larger than expected for a spin-| Dirac particle. 

Earliest definitions of a FF are usually credited to Rosenbluth HRos50l ; in this early reference Rosenbluth 
discussed a model of the proton consisting of a neutron core and a positively charge meson cloud, known then 
as the weak meson coupling model. A high energy electron was expected to penetrate the mesonic cloud and 
to "feel" reduced charges and magnetic moments, e 1 and n'e'. Expressions for such quantities as 
had been derived by Schiff in 1949 HSchi49H . 

In his seminal review paper Hofstadter [Hof56J was the first to relate the results of McAllister and Hofs- 
tadter [McA56J for the ep cross section in elastic scattering at given angle and energy, to the Mott cross section 
for the scattering of a spin | electron by a spin-less proton, ouotu with internal charge density distribution 
p(r), as follows: 

2 



and — 

n e 



;) — ®Mott 



p(r)e i(r -"d 3 r 



volume 



cr M ott\F(q)\ : 



where: 



0" Mott 



2E, 



beam 



cos y 



sin 



4 6» e 



(1) 



(2) 



where E beam and 9 e are the electron incident energy and laboratory scattering angle, respectively, and the 
target mass is infinite. In this early framework a phenomenological FF squared was obtained from absolute 
differential cross section measurements simply as: 
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Figure 1: Fig. 27 in [Hof56], with figure caption 
"The square of the FF plotted against q 2 . q 2 is given in 
units of 10 26 cm -2 . The solid line is calculated for the 
exponential model with rms radii=0.80 x 10~ 13 cm." 



Figure 2: Feynman diagram for the single-photon ex- 
change, or Born term, for elastic ep scattering. 



where q, pbeam and p e are the center-of-mass (CM) momentum transfer, and incident and scattered electron 
momenta, respectively. The historically significant results of these measurements of the proton FF are in 

Fig.m 

2.1.1 The Dirac and Pauli nucleon form factors 

A direct connection between the reduced charge and magnetic moments discussed in [Ros50] and measurable 
observables was first proposed by Clementel and Villi HCle56L who defined FFs on the basis of Rosenbluth's 
discussion of effective charge and magnetic moments, following HSchi49L as F-t (q) = ^ and F 2 (q) = with 
q = 2y 'E beam E e sin y. These FFs were then introduced in experimental papers by Hofstadter and cowork- 
ers HHof56l IMcA56l Hof58], who generalized the "effective" charge and magnetic moment concepts by asso- 
ciating the first with the deviation from a point charge Dirac particle (Dirac FF, Fi), and the second with the 
deviation from a point anomalous magnetic moment (Pauli FF, F 2 ). 

In lowest order, elastic scattering of an electron by the proton is the result of the exchange of a single virtual 
photon of invariant mass squared q 2 = uj 2 — q 2 = —AEb eam E e sin 2 y, (the last step neglects the electron mass), 
where cu = Eb ea m — E e , the energy loss of the electron, and q = pbeam — Pe, the vector momentum change of 
the electron; 6 e is the Lab electron scattering angle. For scattering in the space like region, q 2 is negative. U 

The time-like region, where q 2 is positive, can be accessed for example in e~e + — > pp or pp — > e~e + ; it 
will not be discussed in this review. 

Given the smallness of the fine structure constant a ~ 1/137, it has been common until recently, to neglect 
all higher order terms, except for the next order in a which is treated as a radiative correction, thus implicitly 
assuming that the single photon diagram, corresponding to the Born approximation, is determinant of the 
relation between cross section and FFs; we will revisit this point in section [331 In the single photon-exchange 

'in this review we will use natural units, with energy and mass in GeV, momentum in GeV/c and invariant four-momentum 
transfer squared in (GeV/c) 2 . As is common practice in the literature we will put c=l for convenience and denote momentum 
transfer squared in GeV 2 , although (GeV/c) 2 is understood. 
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process illustrated in Fig. [21 and following the notation of [Rek02|, the amplitude for elastic scattering can be 
written as the product of the four-component leptonic and hadronic currents, i v and J n , respectively: 



iM 



teu(p e )Yu(p beam )} [-ieN(p')T^p / ,p)N(p)] = -t^ 



(4) 



where r 7 ' contains all information of the nucleon structure, u and iV are the electron- and nucleon spinors, 
respectively, g^ u is the metric tensor and k, k', p and p' are the four-momenta of the incident and scattered 
electron and proton, respectively. To ensure relativistic invariance of the amplitude M., can only contain p, 
p' and 7 M , besides scalars, masses and Q 2 . 

As was shown by Foldy HFol52H . the most general form for the hadronic current for a spin |-nucleon with 
internal structure, satisfying relativistic invariance and current conservation is: 
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where Q 2 = —q 2 , is the negative of the square of the invariant mass of the virtual photon in the one-photon- 
exchange approximation in ep scattering, and Fx and F 2 are the two only FFs allowed by relativistic invariance. 
Furthermore, the anomalous part of the magnetic moment for the proton is k p = fi p — 1, and for the neutron 
K n = ^n, in nuclear magneton-units, [i N = ^g, with values k p =1.7928 and n n = —1.9130, respectively; M 
is the nucleon mass. It follows that in the static limit, Q 2 = 0, F lp (0) = 1, F 2p (0) = n p , F ln (0) = 0, -F 2n (0) = 
n n , for the proton and neutron, respectively. 

In the one-photon-exchange approximation Fx(Q 2 ) and F 2 (Q 2 ) are real functions which depend upon Q 2 
only, and are therefore relativistically invariant. When higher order terms with two photons exchange are 
included, there are in general 6 invariant amplitudes, which can be written in terms of 3 complex ones HGuiQ3H . 

The Lab cross section is then: 



da \M\ (E. 
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Following the introduction above, we can now write the standard form for the Lab frame differential cross 
section for ep or en elastic scattering as: 
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where r = Q 2 /4M 2 , and 



2E be 



sin 



T 



1 is the recoil factor. Eq. © is the most 



general form for the cross section, as required by Lorentz invariance, symmetry under space reflection and 
charge conservation. Experimentally, the first separate values for Fx and F 2 were obtained by the intersecting 
ellipse method described by Hofstadter [Hof60J. The early data of Bumiller et al. HBum601 showed that F 2 
decreased with q 2 faster than Fx, even suggesting a diffractive behavior for the proton cross section. Typically 
these results show Fi/F 2 -ratio values which are several times larger than modern values for the proton. 



2.1.2 The electric and magnetic form factors 

Another set of nucleon FFs, F ch and F mag , was first introduced by Yennie, Levy and Ravenhall HYen57h Ernst, 
Sachs and Wali [Ern60J connected F ch and F mag to the charge and current distributions in the nucleon; the 
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interpretation that F ch and F mag measure the interaction with static charge and magnetic fields was given by 
Walecka HWal59H . The following FFs, F ch and F rnag , were defined in HErn60ll : F ch = F 1 - ^jF 2} and F m = 
jjjFx + F 2 . A similar definition of FFs for charge and magnetization, Ge and G M , which is the one in 
use today, was first discussed extensively by Hand, Miller and Wilson [Han 63| who noted that with Ge = 
F\ — tF 2 and G m = F\ + F 2 , the scattering cross section in Eq. © can be written in a much simpler form, 
without interference term, leading to a simple separation method for G Ep and G 2 Mp : 

da ( a V E& ( CQt2( ^ ) \r* i -r 21 ,n rr A cm 
dQ y2E beam sm(f )) E beam y 1 + r J 

Ge p , Gm p , Geu and Gmu are now customarily called the electric- and magnetic Sachs FFs, for the proton and 
neutron, respectively; at Q 2 = they have the static values of the charge and magnetic moments, of the proton 
and neutron, respectively: G Ep (0) = 1, G Mp (0) = H P , G En (0) = and G Mn (0) = ^n- 



2.1.3 Form factors in the Breit frame 

The physical meaning of the electric and magnetic FFs, Ge and G M , is best understood when the hadronic 
current is written in the Breit frame. In that frame the scattered electron transfers momentum qs but no energy 
(lub = 0). Therefore, the proton likewise undergoes only a change of momentum, not of energy, from — 
to +g*s/2; thus Q 2 = q~B 2 . The four components of the hadronic current in the Breit frame are: 

J° = e2Mx n x(F 1 -rF 2 ) = e2M X n xGE, (9) 
J = ie X ^{o x q B )x{Fi + F 2 ) = ie X n (a x q B ) X G M . (10) 

Only in the Breit frame can the electric and magnetic FFs G E and Gm be associated with charge and magnetic 
current density distributions through a Fourier transformation. However the Breit frame is a mathematical con- 
cept without physical reality: there is a Breit frame for every Q 2 value; and above a few GeV 2 , the Breit frame 
moves in the Lab with relativistic velocities, resulting in a non-trivial relation between Breit frame quantities 
and Lab frame quantities: the transformation affects both the kinematics and the structure. A model dependent 
procedure to transform these distributions from the Breit- to the Lab frame has been recently developed by 
Kelly [Kel02J, with interesting results to be discussed later in section H31 



2.2 Rosenbluth form factor separation method 

The Rosenbluth method has been the only technique available to obtain separated values for G\ and G 2 M for 
proton and neutron until the 1990s. The method requires measuring the cross section for eiV scattering at a 
number of electron scattering angles, for a given value of Q 2 ; this is obtained by varying both the beam energy 
and the electron scattering angle over as large a range as experimentally feasible. 

The cross section for ep scattering in Eq. ©, when written in terms of the electric- and magnetic FFs, Ge 
and Gm, takes the following form: 

G 2 m )/(1 + t), (11) 
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and in the notation preferred today, this cross section can be re- written as: 



da ( i 
dn ~ V < 



) * 


G E + — G M 


/(1 + T 


' Mott 







(12) 



where e = [1 + 2(1 + r) tan 2 y] -1 is the virtual photon polarization. 

In early versions of the Rosenbluth separation method for the proton, a correspondingly defined reduced 
cross section was plotted either as a function of cot 2 y HHan63[ IWil64l or cos# e HBer71L For example in 
HHan63L the function R(Q 2 , 9 e ) = [G 2 Ep + rG 2 Mp ] cot 2 f + r(l + r)G 2 Mp was defined. In 1973 Bartel et 
al. chose a form linear in cos 2 y, namely cos 2 y x (^) / {%) Mott HBar73l . Neither of these linearization 
procedures fully disentangles G\ and G 2 Mp . 

The modern version of the Rosenbluth separation technique takes advantage of the linear dependence in e 
of the FFs in the reduced cross section based on Eq. (fT2l) and is defined as follows: 

) 

d ^ J reduced T \ d ^Jexp' J Mott " T 

where (da / dVt) exp is a measured cross section. A fit to several measured reduced cross section values at the 
same Q 2 , but for a range of e-value, gives independently \G\ V as the slope and G 2 M as the intercept, as shown 
in Fig. [3j the data displayed in this figure are taken from [And94 |. 



M e{1+T) (^) /(£) =Gl, + iGl, (13) 
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Figure 3: Demonstration of the Rosenbluth separa- 
tion method based on the data from HAnd94H . The Q 2 
values shown are 2.5 (open triangle), 5.0 (circle) and 
7.0 (filled triangles) GeV 2 . 



Figure 4: Early Rosenbluth separation data for Ge p , 
up to 1973 HBer71l iPrTTTl IBar73l IHan73L but not 
including the 1970 SLAC experiment of Litt et al. 
lEtTOl . 
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2.2.1 Proton form factor measurements 



Figure |4] shows Rosenbluth separation results performed in the 1970's as the ratio Ge p /Ge, where Go is the 
dipole FF given below by Eq. [141 it is noteworthy that these results strongly suggest a decrease of Ge p with 
increasing Q 2 , a fact noted in all four references HBer71[ IPri7 1 1 IBar73[ IHan73L As will be seen in section 
13 .41 the slope of this decrease is about half the one found in recent recoil polarization experiments. Left 
out of this figure are the data of Litt et al. HLit70L the first of a series of SLAC experiments which were 
going to lead to the concept of "scaling" based on Rosenbluth separation results, namely the empirical relation 
I^pGep/Gmp ~ 1. Predictions of the proton FF G Ep made in the same period and shown in Fig. [4] are from 
Refs. HIac73[|Hoh76[|Gar85L all three are based on a dispersion relation description of the FFs, and related to 
the vector meson dominance model (VMD). 
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Figure 5: Data base for Ge p obtained by the Rosen- 
bluth method; the references are I Han63l ILit70[ IPri7Tl 
IBer7 II lBar73l IHan731 lBor75l ISim80l IAnd94l IWal941 
IChr04l[Qat05l. 



Figure 6: Data base for Gm p obtained by the 
Rosenbluth method; the references are BHan631 lJan66 l 
ICow681 iLitTOl iPnTIl IBer711 IHan73l IBar73l IBor75l 
ISil93l lAnd94l IWal94l IChrM [QatOgl . 



A compilation of all Ge p and Gm p data obtained by the the Rosenbluth separation technique is shown in 
Figs. [5] and [6l in these two figures both Ge p and Gm p have been divided by the dipole FF Gd given by: 



G 



D 



l + Q 2 /0.71Ge^ 2 ) 2 



with Ge p — Gd, Gm p — I^ p Gd, and Gmu — I^uGd- 



(14) 



It is apparent from Fig. [5]that the cross section data have lost track of Ge p above Q 2 ~ 1 GeV 2 . It is difficult 
to obtain G\ for large Q 2 values by Rosenbluth separation from ep cross section data for several reasons; first, 
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Q 2 (GeV 2 ) 

Figure 7: The Gm p data follow the pQCD scaling 
law, Gmp ot i, as was first demonstrated in [Arn86|; 
data as in Fig/|6] 



.05 ' 1 1 1 1 1 1 1 1 1 1 1 L - 

0.2 CM 0.5 

tf LfGeWcfj 

Figure 8: Fig. 4 from Ref. HGao03i The 1990 
Platchkov data [Pla90] for Geu, with fits to several 
NN potentials determining the deuteron wave func- 
tion. 



the factor - multiplying G\ in Eq. ([131) automatically reduces the contribution of this term to the cross section 
as Q 2 increases; and second, even at small Q 2 , G 2 M ~ /J%G%, hence the contribution of G\ to the cross section 
is reduced by a factor 7.80. 

In sharp contrast with the situation for Ge p , the Gm p / H p Gd ratios shown in Fig. [6] display excellent internal 
consistency, up to Q 2 = 30 GeV 2 , for the Gj\/ p -values obtained from cross section data; note that the large Q 2 - 
data in HArn86H were obtained without Rosenbluth separation, with the assumption that Ge p = Gm p / 7v the 
ratio Gm p / H p Gd becomes distinctly smaller than 1 above ~ 5 GeV 2 . 

It was first observed by Arnold et al. [Arn86J that the proton magnetic FF, Gm p follows the pQCD predic- 
tion of Brodsky and Farrar [Bro75l, as illustrated in Fig. [71 the pQCD prediction is based on quark counting 
rules. Indeed Q 4 Gm p becomes nearly constant starting at Q 2 = 8 GeV 2 . However, the l/Q 4 -behavior of the 
proton magnetic FF was first mentioned by Coward et al. [Cow68] based on their data extending to 20 GeV 2 ; 
these authors discussed the HQ 4 behavior in light of the vector meson exchange model prevailing at the time 
HSchW67H . 



2.2.2 Neutron electric form factor measurements 

The "neutrality" of the neutron requires the electric FF to be zero at Q 2 = 0, and small at non-zero Q 2 ; 
historically, the fact that the electric FF is non-zero has been explained in terms of a negatively charged pion 
cloud in the neutron, which surrounds a small positive charge HFer47H . 

Early attempts to determine the neutron FF were based on measurements of the elastic ed cross section. 
The scattering by an electron from the spin 1 deuteron requires 3 FFs in the hadronic current operator, for 
the charge, quadrupole and magnetic distributions, Gc, Gq and Gm, respectively. In the original impulse 
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approximation (IA) form of the cross section developed by Gourdin HGou64H , the elastic ed cross section is: 

where A(Q 2 ) = G 2 C (Q 2 ) + lv 2 G 2 Q (Q 2 ) + §7/(1 + V )G 2 M and B(Q 2 ) = §77(1 +77) 2 G^(Q 2 ), with r? = Q 2 /4M^ 
The charge, quadrupole and magnetic FFs can be written in terms of the isoscalar electric and magnetic FFs 
as follows: 

Gc = G e Ce, Gq = G%Cq and Gm = ^{G s M Cs + \G S E C L ), 

where the coefficients Ce, Cq, Cl and Cs are Fourier transforms of specific combinations of the S- and D-state 
deuteron wave functions, u(r) and w(r) HGou64ll . 

The 1971 DESY experiment of Galster et al. HGal71l measured elastic ed cross sections up to 0.6 GeV 2 
with good accuracy and provided a data base for the extraction of Geu', it had been preceded by a series of 
experiments started at the Stanford MARK III accelerator, including Mclntyre and Dhar [McI57J, Friedman, 
Kendall and Gram [Fri60], Drickey and Hand [Dri62J, Benaksas, Drickey and Frerejacque [Ben64J, and Gros- 
setete, Drickey and Lehmann [Gro66a]. On the basis of these data, and using Hamada- Johnston [Ham62J and 
Feshbach-Lomon HFes6711 deuteron wave functions, the following fitting function was proposed in HGal71H : 

G En {Q 2 ) = -7^r^ p (Q 2 ). (16) 
1 + 5. or 

The often quoted Galster fit uses Eq. (fT6l) with Ge v replaced by the dipole FF Go (see Eq. (O). 

The next and last experiment to measure the elastic ed cross section to determine Geu is that of Platchkov 
et al. HPla90l . These data extend to Q 2 of 0.7 GeV 2 , with significantly smaller statistical uncertainties than 
all previous experiments. The data from Platchkov et al. [Pla90| are shown in Fig. [8l The FF A(Q 2 ) is very 
sensitive to the deuteron wave function, and therefore to the NN interaction. Furthermore, the shape of A(Q 2 ) 
cannot be explained by the IA alone. Corrections for meson exchange currents (MEC) and a small contribution 
from relativistic effects were found to significantly improve the agreement between calculations and the shape 
of A(Q 2 ) observed. The authors included the constraint from the slope of the neutron electric FF as determined 
in ne scattering, which at the time was dG En/ dQ 2 =0. 0199 fm 2 [Koe76J; they proposed a modified form of 
the Galster fit using several NN potentials, and including MEC as well as relativistic corrections, of the form 
Geti{Q 2 ) = - a i^G Ep {Q 2 ), corresponding to a slope at Q 2 = 0: dG En /dQ 2 = -afi n /4M 2 . For the Paris 
NN potential for example, a=1.25 and 6=18.3; this fit will be compared with the double-polarization data 
shown later in this review, in Fig. |20l Starting in 1994, all Geu measurements have used either polarization 
transfer or beam-target asymmetry to take advantage of the interference nature of these observables: terms 
proportional to GeuGmu are measured, instead of the G 2 En contribution to the cross section; these experiments 
will be reviewed in section l33l 

2.2.3 Neutron magnetic form factor measurements 

In an early experiment Hughes et al. [Hug65 1 performed a Rosenbluth separation of quasi elastic d(e, e') cross 



sections in the range Q 2 = 0.04 to 1.17 GeV 2 ; they observed non-zero values of G Sn only below 0.2 GeV 2 



2 The isoscalar (Ff) and isovector (F^) Dirac (i = 1) and Pauli (i = 2) FFs are usually defined from the corresponding proton 
and neutron FFs as : Ff = Fi p +Fi n , and F^ — Fi V — Fi n . Analogous relations hold for the Sachs FFs defining G S E = Gep + Geu, 
and G E = Ge p — Geu- 
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but measured G Mn up to 1.17 GeV 2 ; the technique consisted in comparing quasi-elastic ed- with elastic ep 
cross sections. The several experiments following Hughes' can be subdivided into 3 groups: cross section 
measurements in quasi-elastic ed scattering (single arm) HHan73[ lBar731 . which requires large final state in- 
teraction (FSI) corrections at small Q 2 ; elastic ed cross section measurements [Be n64[ lGro66bll : and cross 
section measurements in d(e, e'p)n HBud68[ lDun 66l. or ratio of cross sections d(e, e'ri)p/d(e, e'p)n !lSte66l. 
which is less sensitive to the deuteron wave function, and MEC. All results published prior to 1973 are dis- 
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Figure 9: Early Rosenbluth separation data for Gmu, 
up to 1973 HBen64l |Hug65| IGro66bl ISte66l IDun661 
IBud68llHan73llBar73i 



a 
=1 

\ 
a 
s 



1.25 



1.00 



0.75 



0.50 



0.25 




O RocB2 

• Lun93 
□ Mar93 
B Ank94 
Bru95 
A Ank98 

♦ Kub02 
A Bro05b 



10 



10 



10 



Q 2 (GeV 2 ) 



Figure 10: Recent Gmu data divided by \i n Gr>, 
from cross section data only, starting in 1992 
HRoc92l ILun93l IMar931 IAnk94l lBru95l IAnk98l 
IKub02llBro05ll . 



played in Fig. [9l to be compared with the proton data from the same period in Fig. HI All more recent cross 
section results are in Fig. [10l allowing for a comparison of the progress made in this period for the neutron. 
In Fig. [10] all Gun data obtained from cross section measurements are displayed, including the SLAC ex- 
periments HRoc9 2 , Lun93|, which measured inclusive quasi-elastic ed cross sections. The more recent ELS A 
[Bru95] and MAMI HAnk94[ |Ank98[ [Kub 02 1 experiments are simultaneous measurements of the cross section 
for quasi elastic scattering on the neutron and proton in the deuteron, d(e, e'n)p and d(e, e'p)n; the systemat- 
ics is then dominated by the uncertainty in the neutron detector efficiency; much attention was given to that 
calibration in these experiments. In the ELSA experiment [Bru95| protons and neutrons were detected in the 
same scintillator, and the neutron efficiency was determined in situ with the neutrons from 1 H('j, it). It has 
been argued in [Jou97] that 3-body electro-production contributes significantly and does not necessarily lead 
to a neutron at the 2-body kinematic angle; these data points are shown as <0 in FigJTOl a refutation of these 
arguments is in [Bru97]. In the Mainz experiments the dedicated neutron detector was calibrated in a neutron 
beam at SIN. The new data from Hall B at JLab, HBro05L are shown as filled triangles in Fig. [d for these 
data from Hall B, in addition to the measurements of cross section ratio with the 2 H target, an in-line 1 H target 
was used for an in-situ determination of the neutron counter efficiency via ix + electro-production. 
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Figure 11: Illustration of the quality Figure 12: Illustration of the qual- 
of the dipole fit for Ge p ', data from ity of the dipole fit for Gmp', the data 



Refs. HHan63l iLitTOl IPrTTTl IBer71l 
IHan73l lBar73l IMur74l lBor75l ISim80l 
IWal94llAnd94l in Q 2 range 0.005-2.0 
GeV 2 . 



of Refs. HHan631 lJan66l lLit70l lPrT7Tl 
IBer71l IHan73l IBar73l lBor75l IWal94H 
in the range 0.005-2.0 GeV 2 are in- 
cluded here. 
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Figure 13: Illustration of the qual- 
ity of the dipole fit for Gmu\ the data 
included are the same as in Fig. [TOl 
Refs. HRoc92l lLun93l IMar93l lAnk94l 
IBru95l IAnk98l IKub02l IBro051 in the 
Q 2 -range 0.005-2.0 GeV 2 . 



2.2.4 Rosenbluth results and dipole form factor 

In figures[m[l2]and[l3]the Rosenbluth separation results Ge p , Gm p and Gmu are shown in double logarithmic 
plots for Q 2 < 2 GeV 2 , to emphasize the good agreement of these data with the dipole formula of Eq. O 

Noticeable is the lack of Gm v and Gun data below Q 2 of 0.02 GeV 2 , a consequence of the dominance of 
the electric FF at small Q 2 for the proton, as seen in Eq. (fl~2~l) . 

Although Hofstadter was the first to note that the proton FF data could be fitted by an "exponential model", 
which corresponds to the "dipole model" for FFs in momentum space, it appears that the usage of dividing 
data by Go was introduced first by Goitein et al. HGoi67ll . 

The possible origin of the dipole FF has been discussed in a number of early papers. Within the framework 

VS a v ' s 

of dispersion theory the isovector and isoscalar parts of a FF is written as, G E ' M = Hi 1+Q2 ^ M v,s^ 2 , where 

G V E M , G EM are defined in footnote 2, and M^' s and <\' s are the masses and residua of the isovector-, 
isoscalar vector mesons, respectively. A dipole term occurs when the contribution of two vector mesons with 
opposite residua but similar masses are combined. 
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3 Nucleon form factors from double polarization observables 



It was pointed out in 1968 by Akhiezer and Rekalo [Akh68] that "for large momentum transfers the isolation 
of the charge FF of the proton is difficult" using the elastic ep reaction with an unpolarized electron beam, for 
several reasons: one being G 2 Mp /G 2 Ep > fi 2 at any Q 2 value and the other is that at large Q 2 the contribution 
from the tG 2 m term in G\ + (r/e)G 2 M increases (see Eq.( [T2l) and makes the separation of the charge 
form factor practically impossible. In the same paper the authors also pointed out that the best way to obtain 
the proton charge FF is with polarization experiments, especially by measuring the polarization of the recoil 
proton. Further in a review paper in 1974 Akhiezer and Rekalo [Akh74] discussed specifically the interest of 
measuring an interference term of the form GeGm by measuring the transverse component of the recoiling 
proton polarization in the ep — > ep reaction at large Q 2 , to obtain Ge in the presence of a dominating Gm- hi 
1969, in a review paper Dombey [Dom69| also discussed the virtues of measuring polarization observables in 
elastic and inelastic lepton scattering; however his emphasis was to do these measurements with a polarized 
lepton on a polarized target. Furthermore in 1982 Arnold, Carlson and Gross HArn8 11 emphasized that the best 
way to measure the electric FF of the neutron would be to use the 2 H(e, e'njp reaction. Both a polarized 
target, and a focal plane polarimeter (to measure recoil polarization), have been used to obtain nucleon FFs. 
We discuss below both methods to measure the elastic nucleon FFs, highlighting advantages and disadvantages 
of using polarized target and focal plane polarimeter. 



3.1 Polarization transfer 

Figure [14] shows the kinematical variables for the polarization transfer from a longitudinally polarized electron 
to a struck proton in the one-photon exchange approximation. 




Figure 14: Kinematical variables for polarization 

transfer from a longitudinally polarized electron to a Figure 15: Polarized electron scattering from a polar- 

proton with exchange of a virtual photon. ized target. 

The electron vertex in Fig. [14] can be described by basic Quantum Electrodynamics (QED) rules that 
involves the electron current, £ p , and the proton vertex can be described by QCD and hadron electrodynamics 
involving the hadronic current = x'^F^X- 

For elastic ep scattering with longitudinally polarized electrons, the hadronic tensor, W^ u = J^J*, has four 
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possible terms depending upon the polarization of the target and of the recoil proton: 



= WjfJ + WV(Pi) + WV(P 2 ) + W^(Pi, P a ), (17) 

where the first term in the equation corresponds to unpolarized protons, the second and the third term corre- 
spond to the vector polarization of the initial and the final proton, respectively, and the last term describes the 
reaction when both, the initial and the final protons are polarized. 

Considering the case where only the polarization of the final proton is measured, W^(P) is: 



W^(P) = -TrF^a-P. (18) 



For the scattering of longitudinally polarized electrons off an unpolarized target, in the one-photon-exchange 
approximation, there are only two non-zero polarization components, transverse, P x , and longitudinal, P z \ 
and these components are obtained by calculating the tensors W fll/ (P X ) = ^TrF^F^Ux ~ ItxiGeGm and 
W^{P z ) = \TrF il Fla z ~G\ [ . 

The transformation from Breit to laboratory frame gives following expressions for the polarization compo- 
nents P x and P z in terms of the electric Ge, and magnetic, Gm FF HAkh74[ lArn8 1 1: 

I P X = -2v/r(l + r)G s G M tan| (19) 

I P Z = ^(E beam + E e ) v/r(l + t)G 2 m tan 2 |, (20) 

where E beam and E e are the energy of the incident and scattered electrons, respectively, 9 e is the scattered 
electron angle in the laboratory frame, and I is: 

I = G 2 E {Q 2 ) + -G 2 M {Q 2 ). (21) 
e 

Eqs. (fT9l and (|2Q|) show that the transverse and longitudinal polarization components are proportional to 
G E G M and G\ { , respectively. The ratio G E /G M then can be obtained directly from the ratio of the two 
polarization components P x and P z as follows: 

Ge Px {Ebeam + E e ) 9 e 

= tan — . (22) 

G M P z 2M k 2 j v 

Equation (1221 makes clear that this method offers several experimental advantages over the Rosenbluth sep- 
aration: (1) for a given Q 2 , only a single measurement is necessary, if the polarimeter can measure both 
components at the same time. This greatly reduces the systematic errors associated with angle and beam en- 
ergy change, and (2) the knowledge of the beam polarization and of the analyzing power of the polarimeter is 
not needed to extract the ratio, G E /G M - 



3.2 Asymmetry with polarized targets 

It was pointed out by Dombey [Dom69| that the nucleon FFs can be extracted from the scattering of longitudi- 
nally polarized electrons off a polarized nucleon target. In the one-photon-exchange approximation, following 
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the approach of Donnelly and Raskin HDon86[ |Ras86L the elastic eN (N = p or n) cross section can be 
written as a sum of an unpolarized part and a polarized part; the latter is non-zero only if the electron beam is 
longitudinally polarized: 

a po1 = £ + hA, (23) 

where h is the electron beam helicity, S is the elastic un-polarized cross section given by Eq. (PT21) . and A is 
the polarized part of the cross section with two terms related to the directions of the target polarization. The 
expression for A can be written as [Don86, Ras86J: 



A = -2a Mott tm(e e /2)yJ^-{ y' r[\ - ( I + t) imr(9, /2)\ cos 9" G' u fsiu^ c-os <>G E G } (24) 

where 9* and 0* are the polar and azimuthal laboratory angles of the target polarization vector with q in the u z 
direction and u y normal to the electron scattering plane, as shown in Figure [T51 
The physical asymmetry A is then defined as 

A='-±^ = t (25) 

<7_|_ + <7_ 

where ct + and cr_ are the cross sections for the two beam helicities. 

For a polarized target, the measured asymmetry, A meas , is related to the physical asymmetry, A, by 

Ameas -Pbeam-PtargetA, (26) 

where Pbeam and Pt arg et are electron beam- and target polarization, respectively, and, 



A _Vr(l+T)tan(0 e /2) 



G 2 + iG 2 



sin^* cos (p*G E G M + Jt[1 + (1 + r) tan 2 (0 e /2)] cos^*Gi f 



(27) 



It is evident from Eq. (ITTT ) that to extract G^* the target polarization in the laboratory frame must be perpendic- 
ular with respect to the momentum transfer vector q and within the reaction plane, with 9* = n/2 and (ft* = 0° 
or 180°. For these conditions, the asymmetry A in Eq. (|27l) simplifies to: 



-2y^(lT7ytan(g e /2)g^ 

^iperp — (G B \2 , T ' ^ ' 

^Gm> ^ e 

As (Ge/Gm) 2 is quite small, v4 per . p is approximately proportional to Ge/Gm- In practice, the second term in 
Eq. (|27l) is not strictly zero due to the finite acceptance of the detectors, but these effects are small and depend 
on kinematics only in first order and can be corrected for, so the ratio Ge/Gm is not affected directly. 

The discussion described above is only applicable to a free nucleon; corrections are required if nuclear 
targets, like 2 H or 3 He, are used instead in quasi-elastic scattering to obtain the FFs. 



3.3 Double polarization experiments 

The polarization method, using polarized targets and focal plane polarimeter with longitudinally polarized 
electron beam, has been used to measure both the proton and the neutron e.m. FFs. Below we first describe 
the polarization experiments that measured the proton FFs and next those that measured the neutron FFs. 
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3.3.1 Proton form factor measurements with polarization experiments 



The first experiment to measure the proton polarization observable in ep elastic scattering was done at the Stan- 
ford Linear Accelerator Center (SLAC) by Alguard et al. [ |Alg76 |. They measured the anti-parallel-parallel 



asymmetry in the differential cross sections by scattering longitudinally polarized electrons on polarized pro- 
tons. From their result they concluded that the signs of Ge p and Gm p are the same; they also noted that the 
usefulness of using polarized beam on polarized target is severely limited by low counting rates. 

Next, the recoil polarization method to measure the proton e.m. FF was used at MIT-Bates laboratory 
[Mil98, Bar99]. In this experiment the proton FF ratio Ge p /Gm p was obtained for a free proton and a bound 
proton in a deuterium target at two Q 2 values, 0.38 and 0.5 GeV 2 using polarization transfer from longitudinally 
polarized electron to the proton in the target, and measuring the polarization of the recoiling proton with a focal 
plane polarimeter (FPP). The conclusion from these measurements was that the polarization transfer technique 
showed great promise for future measurements of Ge p and Geu at higher Q 2 values. 

The ratio Ge p /Gm p in elastic 1 H(e,e'p) was also measured at the MAMI in a dedicated experiment 
UPosOlL and as a calibration measurement [DieOl] at a Q 2 of ^0.4 Gev 2 . The ratio values found were in 
agreement with other polarization measurements as well as Rosenbluth measurements. Most recently the 
BLAST group at MIT-Bates HCra06H measured the ratio G Ep /G Mp at Q 2 values of 0.2 to 0.6 GeV 2 with high 
precision. 

Starting in late 1990's, the proton FF ratios G Ep /Gm p were measured in two successive dedicated exper- 
iments in Hall A at JLab for Q 2 from 0.5 to 5.6 GeV 2 HJonOOl |Gay02[ IPun 051. Other measurements were 



also conducted in Hall A [GayOTj |Str03[ |Hu06H at lower Q 2 values, as calibration measurements for other 



polarization experiments, and one measurement in Hall C [Mac L06ll . 

In the first JLab experiment the ratio, G Ep /Gm p , was measured up to Q 2 of 3.5 GeV 2 [JonOOl |Pun05H . 
Protons and electrons were detected in coincidence in the two high-resolution spectrometers (HRS) of Hall A 
[Alc04J. The polarization of the recoiling proton was obtained from the asymmetry of the azimuthal distribu- 
tion after the proton re-scattered in a graphite analyzer. 

The ratio, Ge p /Gm p , was measured at Q 2 = 4.0, 4.8 and 5.6 GeV 2 with an overlap point at Q 2 = 3.5 GeV 2 
[Gay02[ , in the second JLab experiment. To extend the measurement to higher Q 2 , two changes were made 
from the first experiment. First, to increase the figure-of-merit (FOM) of the FPP, a CH 2 analyzer was used 
instead of graphite; hydrogen has much higher analyzing power [Spi83, Mil77] than carbon [Che95]. Second, 



the thickness of the analyzer was increased from 50 cm of graphite to 100 cm of CH 2 to increase the fraction of 
events with a second scattering in the analyzer. Third, the electrons were detected in a lead-glass calorimeter 
with a large frontal area, to achieve complete solid angle matching with the HRS detecting the proton. At the 
largest Q 2 of 5.6 GeV 2 the solid angle of the calorimeter was 6 times that of the HRS. 

Proton polarimeters are based on nuclear scattering from an analyzer material like graphite or CH 2 ; the 
proton-nucleus spin-orbit interaction and proton-proton spin dependent interaction results in an azimuthal 
asymmetry in the scattering distribution which can be analyzed to obtain the proton polarization. The detection 
probability for a proton scattered by the analyzer with polar angle $ and azimuthal angle (p is given by: 

<p) = € -^- (1 ± AyiPtewup- P/ pp cos^)) , (29) 

where ± refers to the sign of the beam helicity, Pf pp and Pf pp are transverse and normal polarization compo- 
nents in the reaction plane at the analyzer, respectively, and e(i9, ip) is an instrumental asymmetry that describes 
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a non-uniform detector response. Physical asymmetries are obtained from the difference distribution of f ± , 



A = {ft - /D/2 = ^ (A^ sm ^ - AyPte cos ^) , (30) 
and the sum distribution of f ± separates the instrumental asymmetries e$, 

Ei = (/+ + /r)/2 = ^- (3D 

The values of the two asymmetries at the FPP, A y Pf pp and A y Pf pp , can be obtained by Fourier analysis of the 
difference distribution Di, however to calculate the ratio Ge p /Gm p , the proton polarization components P x 
and P z are needed at the target. 

As the proton travels from the target to the focal plane through the magnetic elements of the HRS, its spin 
precesses, and therefore the polarization components at the FPP and at the target are different. The hadron 
HRS in Hall A consists of three quadrupoles and one dipole with shaped entrance and exit edges, as well as a 
radial field gradient. The polarization vectors at the polarimeter, pf pp , are related to those at the target, hP, 
through a 3-dimensional spin rotation matrix, (S), as follows: 

pfpp 
y 

pfpp 

X 

pfpp 

z 

The spin transport matrix elements SV, can be calculated using a model of the HRS with quadrupoles, fringe 
fields, and radial field gradient in the dipole, for each tuning of the spectrometer setting, and event by event 
with the differential-algebra-based transport code COSY [Ber95 1. The spin transport method to obtain the two 
asymmetries at the target, hA y P x and hA y P z , was developed by Pentchev and described in detail in [Pe n03L 
and discussed in [Pun05J. The ratio Ge p /Gm p was calculated from the two asymmetries at the target from 
Eq. (|22)) . The fact that both beam polarization, and polarimeter analyzing power cancel out of this equation 
contributes to the reduction of the systematic uncertainties, however their values do influence the statistical 
errors. 

The most recent acquisition of the FF ratio Ge p /Gm p has been made at a Q 2 of 1.51 GeV 2 by measuring 
the beam-target asymmetry in an experiment in Hall C at JLab in elastic ep scattering [Jon06J. This is the 
highest Q 2 at which the Ge p /Gm p ratio has been obtained from a beam-target asymmetry measurement. 
The results from the two JLab experiments [ Gay02[|Pun05i and other polarization measurements [Mil98, 




Gay01[ IPosOll IDieOll IStr03l lHu06l MacL06, Jon06J, are plotted in Fig. [16] as the ratio fi p G Ep /G Mp versus 



Q 2 . All data show only the statistical uncertainty; the systematic uncertainty for the data of [GayOfj IPun05l 
are shown separately as a polygon; they are typical for all polarization data obtained in Hall A at JLab. The 
new asymmetry data from BATES [Cra06[ are not in this figure as they are in the range of Q 2 - values smaller 
than 0.6 GeV 2 ; they appear in Fig. [221 As can be seen from figure [161 data from different experiments are in 
excellent agreement and the statistical uncertainty is small for all data points; this is unlike Ge p obtained from 
cross section data and shown in Fig. [51 where we see a large scatter in results from different experiments as 
well as large statistical uncertainty at higher Q 2 values, underlining the difficulties in obtaining Ge p by the 
Rosenbluth separation method. 

The results from the two JLab experiments [JonOO, Pun031 pay02[ showed conclusively for the first time 



a clear deviation of the proton FF ratio from unity, starting at Q 2 ~ 1 GeV 2 ; older data from [Ber7TJ IPri7 1 1 
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Bar73 , Han73J showed such a decreasing ratio, but with much larger statistical and systematic uncertainties, 
as seen in Fig. HI The most important feature of the JLab data, is the sharp decrease of the ratio [i p Ge p /Gm p 
from 1 starting at Q 2 m 1 GeV 2 to a value of ~ 0.28 at Q 2 = 5.6 GeV 2 , which indicates that Ge p falls faster 
with increasing Q 2 than Gm p - This was the first definite experimental indication that the Q 2 dependence of 
G Ep and G Mp is different. If the [i p G Ep l GMp-ratio continues the observed linear decrease with the same slope, 
it will cross zero at Q 2 ~ 7.5 GeV 2 and become negative. 

In Fig. [T7]all the ratio data obtained from the Rosenbluth separation method are plotted together with the re- 
sults of the two JLab polarization experiments. There are recent proton FF results obtained with the Rosenbluth 
separation method from two JLab experiments [Chr04, Qat05]; these results agree with previous Rosenbluth 
results HLit70[ lBer7 1 [ IPriTTl |Wal94l I And94| and confirm the discrepancy between the ratios obtained with the 
Rosenbluth separation method and the recoil polarization method. The two methods give definitively different 
results; the difference cannot be bridged by either simple re-normalization of the Rosenbluth data HArr03H . or 
by variation of the polarization data within the quoted statistical and systematic uncertainties. This discrep- 
ancy has been known for the past several years and is currently the subject of intense discussion. A possible 
explanation is the hard two-photon exchange process, which affects both cross section and polarization transfer 
components at the level of only a few percents; however, in some calculations [Afa0fl lBlu03ll the contribution 
of the two-photon process has drastic effect on the Rosenbluth separation results, whereas in others it does not 
[Bys06|; in either case it modifies the ratio obtained with the polarization method by a few percent only (this 
will be discussed in section [331) . There are several experiments planned at JLab [Sul04, Arr05J to investigate 
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the two-photon effects in the near future. 



3.3.2 Neutron electric form factor measurements with polarization experiments 

Measurements of the FFs of the neutron are far more difficult than for the proton, mainly because there are 
no free neutron targets. Neutron FF measurements were started at about the same time as for the proton, but 
the data are generally not of the same quality as for the proton, especially in the case of the electric FF of 
the neutron; the Q 2 range is limited also. The early measurements of the FFs of the neutron are discussed 
in sections [2.2.21 and 12.2.31 in this section we discuss only measurements made with longitudinally polarized 
electron beams on polarized 2 H- or 3 He-targets, and polarization transfer in the 2 H(e, e'n)p reaction. We start 
with the measurements of the charge FF, Geh, and proceed then to the relatively few measurements of the 
magnetic FF, G Mn - 

The first measurement of the charge FF of the neutron, Ge u , by the polarization method was made at 
MIT-Bates using the exclusive 2 H(e, e'n)p reaction [Ede94j. The advantage of polarization measurements on 
the deuteron in the quasi free kinematics is that the extracted neutron FF is quite insensitive to the choice of 
deuteron wave functions, and also to higher order effects like final state interaction (FSI), meson exchange 
currents (MEC) and isobar configurations (IC), when the momentum of the knocked out neutron is in the 
direction of three-momentum transfer q HAre87[ |Rek89, Lag91 1. 



For a free neutron the polarization transfer coefficient P x is given by Eq. (fT9l) . The relation between polar- 
ization transfer coefficient P x , the beam polarization, P e , and the measured neutron polarization component, 
P x , is P' x = P e P x - The FF G En was extracted at a Q 2 of 0.255 GeV 2 in this experiment from the measured 
transverse polarization component P' x of the recoiling neutron, and known beam polarization, P e . This early 
experiment demonstrated the feasibility of extracting Geu from the quasi-elastic 2 H(e, e'n)p reaction with the 
recoil polarization technique, with the possibility of extension to larger Q 2 values. 

Next, this same reaction 2 H(e, e'n)p was used to determine Geu at MAMI [ Her99[ IOst99l by measuring 
the neutron recoil polarization ratio P x /P z , at a Q 2 of 0.15 and 0.34 GeV 2 . The ratio P x /P z is related to 
Ge/Gm as shown in Eq. (|22|) . The measurement of the ratio, P x /P z , has some advantage, as discussed 
earlier for the proton, over the measurement of P x only, because in the ratio the electron beam polarization 
and the polarimeter analyzing power cancel; as a result the systematic uncertainty is small. In yet another 
experiment at MAMI the ratio of polarization transfer components, P x /P z , was measured using the same 
reaction 2 H(e, e'n)p and the electric FF G En was obtained at Q 2 = 0.3, 0.6 and 0.8 GeV 2 HGla05ll ; Glazier et 
al. concluded that their results were in good agreement with all other Geh double-polarization measurements. 

The experiment at JLab by Madey et al. HMad03[ IPlaQ5B obtained the neutron FF ratios Ge/Gm at Q 2 
values of 0.45, 1.13 and 1.45 GeV 2 using the same method of measuring the recoil neutron polarization com- 
ponents P x and P z simultaneously, using a dipole with vertical B-field to precess the neutron polarization in 
the reaction plane, hence obtaining directly the ratio GeuIGmu- The best-fit values of Gm u were used to 
calculate values of Ge u from the ratio measurements. This is the first experiment that determined the value of 
Geu with small statistical and systematic uncertainty and at the relatively high Q 2 values up to 1.45 GeV 2 . 

Passchier et al. [Pas99| reported the first measurement of spin-correlation parameters A^ d at a Q 2 of 
0.21 GeV 2 in 2 H (e, e'n)p reaction at NIKHEF; this experiment used a stored polarized electron beam and an 
internal vector polarized deuterium gas target; they extracted the value of Ge u from the measured sideways 
spin-correlation parameter in quasi-free scattering. 

Experiment E93-026 at JLab extracted the neutron electric FF at Q 2 = 0.5 and 1.0 GeV 2 HZhuOHlWar04ll 
from measurements of the beam-target asymmetry using the 2 H(e, e'n)p reaction in quasi elastic kinematics; 
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in this experiment the polarized electrons were scattered off a polarized deuterated ammonia (ND 3 ) target. 
This experiment was the first to obtain Geu a t a relatively large Q 2 using a polarized target. 

Blankleider and Woloshyn, in a paper in 1984 HBla84L proposed that a polarized 3 He target could be used 
to measure G En or Gyin- They argued that the 3 He ground state is dominated by the spatially symmetric S-state 
in which the two proton spins point in opposite directions, hence the spin of the nucleus is largely carried by 
the neutron. Therefore, the 3 He target effectively serves as a polarized neutron target; and in the quasi-elastic 
scattering region the spin-dependent properties are dominated by the neutron in the 3 He target. 

There were experiments in the early 1990's at MIT-Bates Laboratory that used a polarized 3 He target and 
measured the asymmetry with polarized electrons in spin-dependent quasi-elastic scattering BJon91[ |Tho92H , 
and extracted the value of Geu using the prescription of Blankleider and Woloshyn [Bla84J, at a Q 2 =0.\6 
and 0.2 GeV 2 . However, Thompson et al. [Tho92] pointed out that significant corrections are necessary at 
Q 2 =0.2 GeV 2 for spin-dependent quasi elastic scattering on polarized 3 He according to the calculation of Laget 
[Lag91 ]; hence no useful information on Geu could be extracted from these measurements; but Thompson et 
al. concluded that at higher Q 2 values the relative contribution of the polarized protons becomes significantly 
less and a precise measurements of Geu using polarized 3 He targets will become possible. 

Starting in the early 1990's, the neutron electric FF G En has been obtained in several experiments at MAMI, 
by measuring the beam-target asymmetry in the exclusive quasi-elastic scattering of electrons from polarized 
3 He in the 3 He(e, e'n)pp reaction [Mey94, Bec99l |Roh99l |Ber03H . All the Ge u data from polarization exper- 
iments are shown in Fig. In the first of these experiments, at MAMI, [Mey94|, Geu was obtained at Q 2 
= 0.31 GeV 2 . In following experiments at MAMI, G En was extracted at Q 2 of 0.35 GeV 2 HBec991 and 0.67 
GeV 2 HRoh991lBer03l using the same reaction. The 0.35 GeV 2 point of HBec99D was later corrected in HGolOll . 
based on Faddeev solutions and with some MEC corrections. The large effect of these corrections is illustrated 
in Fig. [T8]with the dashed line connecting the open diamonds. The size of these corrections is expected to 
decrease with Q 2 , although the corrections become increasingly difficult to calculate with increasing Q 2 s. 



3.3.3 Neutron magnetic form factor measurements with polarization experiments 

Only two experiments have obtained the magnetic FF of the neutron, Gmu, from polarization observables; both 
experiments used a polarized 3 He target. The first experiment at the MIT-Bates laboratory, extracted Gmu from 
the measured beam-target asymmetry in inclusive quasi-elastic scattering of polarized electrons from polarized 
3 He target at Q 2 of 0.19 GeV 2 HGao94M ; the uncertainty on Gun was dominated by the statistics, with a 
relatively small contribution from model dependence of the analysis. The second JLab experiment obtained 
Gmu f° r Q 2 values between 0.1 and 0.6 GeV 2 , by measuring the transverse asymmetry in the 3 He(e,e') 
reaction in quasi-free kinematics HXuOOl IXu0 3 , And07]. The values of Gun were extracted in the plane wave 
impulse approximation (PWIA) at Q 2 of 0.3 to 0.6 GeV 2 , and from a full Faddeev calculation at Q 2 of 0.1 and 
0.2 GeV 2 . The authors of this paper asserted that the PWIA extraction of Gmu is reasonably reliable in the 
Q 2 range of 0.3 to 0.6 GeV 2 ; however, a more precise extraction of G Mn requires fully relativistic three-body 
calculations. The Gm u values from both experiments are shown in Fig. [191 



3.4 Discussion of the form factor data 

Probably the most important advance in the characterization of the FFs of the nucleon made in the last 10 years 
has been the realization that the so-called "scaling"-behavior of the proton FFs: 

Gep~Gm p IHp ~ Gd, (32) 
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Figure 19: The data for Gmu from the two polariza- 
tion experiments which have been performed so far, at 
Bates MGao941 and JLab IAnd071 . respectively. 



was limited to values of Q 2 smaller than 2 GeV 2 . The recoil polarization data obtained at JLab in 1998 and 
2000 proved beyond any doubt that for Q 2 -values larger than 2 GeV 2 , Ge p decreases faster than Gm p / Hp with 
a slope of -0.14 per GeV 2 . What we now have are distinctly different Q 2 -dependences for Ge p and Gm p ', in 
the Q 2 region investigated so far, the scaling-behavior is violated by a factor of 3.66 to']® at 5.54 GeV 2 . The 
deviation of Ge p from the dipole FF is illustrated in Fig. [22l where only polarization results are shown. Of 
course it was well known that the dipole FFs, when Fourier transformed, produce unphysical distributions of 
charge or magnetization, with a discontinuity at zero radius. Nevertheless there were valid reasons, to believe 
that the dipole FF discussed in 12.2.41 may actually describe the FFs Ge p , Gm p and Gmb of the nucleon. The 
data no longer support such expectations, as can be concluded by comparing the results in Figs. I5landl22l 

The discrepancy is related to the techniques used: all Rosenbluth separation of cross section data including 
the 2 new measurements from JLab [Chr04, Q at05[ give n P G Ep /G Mp ratios close to the scaling behavior, 



except the early data shown in Fig. 3; all recoil polarization results for the same ratio are clustered along an 
approximately straight line versus Q 2 , with a best fit valid above Q 2 ~ 0.4 GeV 2 given by: 

H P G Ep /G Mp = 1.0587 - 0.14265 Q 2 . (33) 

A number of observations relative to this difference in results follows. First, there is one well established dif- 
ference between the two techniques, cross section versus recoil polarization, and it is the relative importance of 
the radiative corrections required for them, as discussed in detail in section [331 The total radiative corrections 
as routinely calculated in cross section measurements is typically 10 to 30%, and the corrections are strongly e 
dependent; this e dependence affects primarily the results for G 2 Ep , and for increasing Q 2 the accuracy require- 
ment for the correction becomes very demanding. Second, polarization observables, in recoil polarization 
or target asymmetry measurements, being ratios of cross sections, are only minimally affected by radiative 
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corrections, and the ratio Ge/Gm even less being a ratio of ratios. Nevertheless polarization data ultimately 
will require radiative corrections, particularly as experiments continue into the domain of yet larger Q 2 . So is 
the discrepancy between Rosenbluth and polarization data entirely due to inaccuracy or incompleteness in the 
radiative correction? An immediate consequence of the previous statements is that radiative corrections for 
elastic ep scattering in general have to be reexamined, as in their presently practiced form they are unable to 
reconcile the cross section results with polarization results. 

Encouraging progress has been made including the one process certainly neglected in all previous ra- 
diative corrections, the exchange of two photons, neither one of them "soft" (this will be further discussed 
in section [331 ) . Several calculations [Gui03, Afa05al |Blu03ll suggest that this one diagram may contribute 
significantly to the e-dependence of the cross section; other considerations lead to the conclusion that the con- 
tribution from the two-photon term is too small at the Q 2 -values of interest JBys06[ , and/or leads to a definite 
non-linearity in the Rosenbluth plot which has not been seen in the data so far HTom05ll . 

Following the publication of the JLab recoil polarization Ge p / Gm p ratios up to 5.54 GeV 2 , the entire cross 
section data base for the proton has been reanalyzed by Brash et al. HBra02ll . leaving all data above Q 2 = 1 
GeV 2 out, using the data from [JonOO] Gay02 1 above this value of Q 2 , and allowing for relative renormalization 
of all cross section data so as to minimize the x 2 of a global fit for Gm p - The fitting function is the inverse of a 
polynomial of order 5. The renormalized values of Gm p show less scatter than the original data base, and the 
net effect of imposing the recoil polarization results is to re-normalize all Gm p data upward by 1.5-3% when 
compared with the older Bosted parametrization [Bos95 |, as shown in Fig. |23] 

Another useful fit to the nucleon FFs which gives a good representation of the data is the one by Kelly 
[Kel04j. This fit uses ratios of polynomials with maximum powers chosen such that Ge p , Gm p and Gmu 
have the asymptotic l/Q 4 behavior required by pQCD; in [Kel04j Geu was also re-fitted with a Galster FF, as 
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Figure 22: Polarization data presented as Ge p /Gd, 
where Ge p is obtained from the ratio Ge p /Gmp 
obtained from polarization data in MPun05[ Gay 02 [ 
ICra06i multiplied by G Mp from the Kelly fit HKel04i 
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shown in Fig.l20l 

In Figs. [20] and [22 we compare all the data available for Geu and Gmu, obtained from cross section and 
polarization observables. The Geu data obtained in double polarization show reasonable consistency above 0.5 
GeV 2 ; they are systematically higher than the older cross section results shown in Fig. [20]by the 3 Platchkov 
fits HPla901 . The revision by Kelly [Kel04] of the Galster fit HGal71ll gives an excellent representation of the 
data available today. 

Recently, more and better data have been obtained for Gs n , exclusively by the polarization method, either 
recoil polarization transfer or target asymmetry, with deuterium and 3 He targets and up to Q 2 =\.5 GeV 2 . 
No drastic change of the general behavior of Geu has been observed to this point in time. There is a new 
measurement of Geu at JLab up to a Q 2 of 3.4 GeV 2 HCat031 . but the data have yet to be analyzed. In general 
all polarization data for Geu have given results larger than those obtained from elastic scattering; these earlier 
data required considerable nuclear structure corrections, as illustrated in Fig. [TH the sensitivity to the deuteron 
wave function, therefore to the NN potential used, was extensively discussed at the time in [Pl a90l . 

The data for G Mn come mostly from cross section measurements, except two polarization measurements, 
using polarized 3 He target, one at MIT-Bates for low Q 2 with large uncertainty [Gao94| and the other a recent 
measurement at JLab HXuOOI IXu03l |And07i The most recent Hall B results HBro05L which extend to Q 2 of 
nearly 5 GeV 2 , and used quasi-elastic scattering on deuterium, reveal some internal inconsistency in the data 
base near 1 GeV 2 as shown in Fig. [2H as can be seen in Fig. [2T] there is some disagreement between the 
results of different experiments in the Q 2 range of 0.3 to 1.5 GeV 2 . . These measurements will be extended to 
14 GeV 2 after the JLab upgrade to 12 GeV; similarly, the measurement of Ge p /Gm p will be continued to 13 
GeV 2 after the upgrade. 

Several experiments are planned at JLab to resolve the dichotomy in the Ge p /Gm p ratio. One experi- 
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ment will measure the ratio of the e~p and e + p cross sections, which determines directly the real part of the 
two-photon amplitude [Afa04j. Another experiment will measure the ratio Ge p /Gm p at fixed Q 2 =2.5 GeV 2 
HSul04ll . as a function of e, to detect the two hard photon contribution as a variation of this ratio; non constancy 
would be related to the real part of the two-photon amplitude. A third experiment will be a high statistics 
search of non-linearity in the Rosenbluth plot in ep-scattering, which should also reveal the contribution of 
the two-photon process ( HArr Q51). Measurements of the induced polarization in ep — > ep (a byproduct of the 
experiment from HSul04|Q , and of the single spin target asymmetry in quasi elastic scattering on the neutron 
in 3 He^ (e, e') for target polarization normal to the reaction plane [Ave05], will measure the imaginary part of 
the two-gamma contribution. The transverse beam spin asymmetry in ep has been measured at Bates [WelOlJ 
and MAMI [Maa05J; it too originates from the imaginary part of the two-photon contribution. 

3.5 Rosenbluth results and radiative corrections 

All cross section measurements have been single arm experiments, (e, e'), except three early experiments at 
Cambridge IPriTTl IHan73ll and DESY HBar73H in which both proton and electron were detected, (e, e'p), and 
the most recent one of Qattan et al. [ |Qat05| , in which only the proton was detected, (e,p). In all cases, 
measured raw cross sections need to be corrected for QED processes to first order in a ~ 1/137, before 
accessing the cross section corresponding to one-photon-exchange, or Born term. Only to the extent that these 
corrections remain relatively small, can one hope to obtain the Born term FFs G\ and G 2 M , which are functions 
of Q 2 only, using the Rosenbluth method. 

The effect of the radiative correction on the cross section is typically in the range 10-30%; what is important 
however, is the fact that overall, the radiative corrections are e-dependent; i.e. they affect the slope of the 
Rosenbluth plot. Although radiative corrections have been applied to all data taken after 1966 using the 
"recipe" of Tsai [Tsa61J, Mo and Tsai [MoT69| and HTsa71L not all corrections were applied in all data sets. 
This point was recently reviewed by Arrington [Arr03J, who reanalyzed some of the cross section data; the fit to 
the re-analyzed data is included in Fig. [FT] Furthermore, in the references HTsa61[|Mo T69, Tsa71J the effect of 
the structure of the nucleon was ignored, and a number of approximations were made. In more recent work on 
radiative corrections, Maximon and Tjon [MaxOO] have included the structure of the proton by introducing the 
proton FF, and they also eliminated some of the soft-photon approximations made by HTsa61[|Mo T69, Ts a71H . 
In the current energy range of JLab, the difference for 5, the radiative correction, used in ^ = (1 + 5) dc r d ^" 
up to corrections of order a 3 , between the older and the new calculation is at the level of several %. 

The various internal radiative correction diagrams involving the electron are shown in Fig. [211 The first 
order virtual radiative processes are the vertex diagram b), the photon self-energy diagram c) and the two self- 
energy diagrams for the electron d); the first order real radiative processes include emission of a real photon 
by either the initial or the final electron diagram e). Similarly diagrams for the proton include bremsstrahlung 
a), vertex b) and proton self energy c), shown in Fig. [25l Two-photon exchange is shown as diagram e). 
In addition there are external radiative corrections due to the emission of real photons by the incoming and 
scattered electrons in the material of the target, as well as energy loss by ionization. 

The virtual part of the internal radiative corrections depend exclusively upon Q 2 , thus it generates no 
e-dependence, hence does not modify the value of G\ , but modifies the value of G 2 Mp directly. 

The radiative correction for real photon emission (bremsstrahlung) is energy, and therefore e dependent, 
and it also results in a changed value of Q 2 . In general the scattered electron energy spectrum is integrated 
up to a maximum energy loss which is kept below the pion threshold. The correction is different for different 
experiments; it depends on the procedure used to integrate over the scattered electron energy, or missing mass 
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Figure 24: Born term and lowest order radiative 
correction graphs for the electron in elastic ep. 



Figure 25: Lowest order radiative correction for 
the proton side in elastic ep scattering. 



squared spectrum. 

The contributions due to real photon emission by the initial and final proton, as well as the proton vertex 
and two-photon exchange with one soft- and one hard photon are relatively small, but strongly e-dependent. 

The external part of the radiative corrections includes only real photon emission by the incident and scat- 
tered electron, and is not coherent with the ep interaction. Although the correction for the incoming electron 
in the target is energy independent, and it can be averaged to a value at the center of the active area of the 
target for all kinematics of a given experiment, the correction for the scattered electron in the target depends 
directly upon the target length and diameter which determines the amount of target material traversed, and 
therefore the scattering angle. As the desired range of e values is obtained by changing the electron scattering 
angle, this correction has e-dependence. For the data of Andivahis et al. [And94j the external corrections are 
one fourth to one half as large as the internal corrections from the smallest to the largest e- values as shown in 
Fig. [261 The calculation of the external correction requires information on the spectrometer acceptance and 
on the target geometry, and is an integral part of the analysis of the data; it cannot be repeated on the basis of 
published data. However, it is potentially a significant source of uncertainty in the e-dependence of the total 
radiative correction. 

To gain some appreciation of what term might be most strongly affecting the final result of the radiative 
correction, we show the values 5 rea i + 5 virtua i, graphs b), c), d) and e) in Fig. [2U 5 12 from graphs a), b),c) and 



d) in Fig. [25j and 8 e 



rtemai f° r the condition of the Andivahis experiment [And94] separately in Fig [261 The 
curve in Fig. [26] labeled 5 int+ ext determines the overall correction. Its slope versus e is due to the combined 
effect of the real and external-contributions, with the proton contribution reducing it somewhat; 5 V i rtU ai has no 
e-dependence. 

The importance of calculating the contributions to the radiative correction which are e-dependent accurately 
is illustrated in Fig. [27j Shown in this figure are reduced cross sections defined in terms of the (^) re duced from 
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Figure 28: The Q 2 = 5 GeV 2 data of 
MAnd94|| , with their best fit (thick solid 
line), compared with the results of the 
various calculations of the radiative 
correction, taken from Refs. [VdhOO, 
IAfa05al |Bys06[ . Also shown is the 
slope from the JLab polarization data. 



Eq. (H3) as: 



da\ .^,2 _ 6 G\ v ^ Gj 

d^ 1 ) reduced T G D G D 



^=(^1 + (34) 



as a function of e for the data of Andivahis et al. [And94|. If both FFs are functions of Q 2 only, the intercept of 

a straight line fit is G\ Ip /G 2 D ~ fi 2 , and the slope is Most noticeable in this figure is the negative slope 

of the uncorrected data, above Q 2 =3 GeV 2 . This figure dramatically illustrates the importance of the radiative 
correction and gives a measure of the accuracy that is required to obtain the FFs with the desired accuracy. 
The final value of G\ obtained from cross section data depends directly upon the value and the accuracy of 
the e-dependent part of the radiative correction. Note that the radiative corrections for the data of [And94| 
were made following Mo and Tsai [Tsa61, MoT69l lTsa71L with the additional corrections introduced in Ref. 
HWal94I 

More recently Maximon and Tjon [MaxOOJ have reconsidered the radiative correction calculation, and 
included additional terms with explicit emphasis of the hadronic effects. A similar reexamination of the Mo- 
Tsai procedure was made by Vanderhaeghen et al. [VdhOO] in the process of a detailed calculation of radiative 
corrections for virtual Compton (VCS). Also recently Ent et al. HEntOlU and Weissbach et al. [Wei04] have 
published improvements and detailed studies of the radiative correction calculation technique for coincidence 
experiments (e,e'p). 

Most recently Bystritskiy et al. [ Bys06[ have calculated the radiative corrections for elastic ep scattering 
using the Drell-Yan electron structure function approach; no co-linearity approximation is made in such a cal- 
culation, but the proton vertex corrections have not been included so far; the diagram with two hard photons has 
been approximatively calculated using both nucleon and A intermediate states and was found to make a negli- 
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gible contribution. The results of [Bys06| suggest that hard bremsstrahlung may cause the difference between 
the Rosenbluth and polarization results. Usual bremsstrahlung calculations are for soft bremsstrahlung, where 
the emitted photon energy is kept only to linear order in denominators and entirely omitted in numerators. Soft 
bremsstrahlung multiplies all amplitudes by the same factor and does not, for a relevant example, change the 
slope on a Rosenbluth plot. If one makes no approximations in the photon energy, there can be different effects 
on different spin amplitudes. Thus the claim is that emitted photons that are energetic enough to affect the spin 
structure of the calculation but still small enough to escape detection, give rise to the difference between the 
two methods of measuring Ge p /Gm p - A contrasting numerical claim is that hard bremsstrahlung effects are 
noticeable and helpful in reconciling the Rosenbluth and polarization experiments, but are not decisive, see 
Ref. [Afa05bJ. These contrasting claims clearly need to be sorted out, but an independent reexamination is not 
available as of this writing. 

The effect of these new radiative corrections is illustrated in Fig. [28] The dashed-dot line is obtained from 
the Q 2 =5 GeV 2 , uncorrected data from HAnd94| , applying the radiative correction calculated with the code 
of HVdhOOL with the same energy cuts as used in the original data. This correction is 2.5% smaller than the 
one in [And94l at e=l. The two-photon calculation result shown is obtained by removing the soft part of the 
two-photon contribution, and replacing it by the GPD based calculation of [Afa05a]; the result is then refitted 
with a straight line (dotted line). The resulting value at e=l is 4.5% smaller than the original correction. The 
long dashed line represents the results of [Bys06], after correcting the experimental data points HTom061 . us- 
ing the same energy cuts as in the data, and refitting with a straight line (dashed line); these results are almost 
identical to the ones obtained with two-photon correction [Afa05a|. The slope calculated from the fit to the 
JLab recoil polarization data, for Q 2 =5 GeV 2 , from Eq. (|33l ), is shown as the short dashed line. The value at 
e=l is 6% smaller than that of the original Rosenbluth data of Ref.[And94|. In Fig. [28] all fits are drawn with 
a renormalized value of cr^ at e=0, to emphasize the differences in slope, which determine Ge p ', based on the 
recoil polarization results the contribution of Ge p to the cross section at Q 2 =5 GeV 2 is ~1%. 

All three corrections are different, and each one of them brings the Rosenbluth results closer to the recoil 
polarization results, indicating that present uncertainties in the calculations of the radiative corrections of the 
cross section are at the level of several %s. 

In HArr04ll . the effect of the Coulomb distortion of the incoming and outgoing electron waves on the ex- 
traction of the proton FFs was studied. Coulomb distortion corresponds to the exchange of one hard and one 
(or several) soft photons. It was found that it does yield an e dependent correction to the elastic electron-proton 
cross sections. Although it reduces the cross sections, its magnitude is too small to explain the discrepancy be- 
tween Rosenbluth and polarization methods. It is however straightforward to calculate and should be included 
in the data analysis. 

Following the important discrepancy between the determinations of Ge p /Gm p using the polarization trans- 
fer and Rosenbluth techniques, the role of two hard photon exchange effects, beyond those which have already 
been accounted for in the standard treatment of radiative corrections has been studied. A general study of 
two- (and multi)-photon exchange contributions to the elastic electron-proton scattering observables was given 
in HGui031l . In that work, it was noted that the interference of the two-photon exchange amplitude with the one- 
photon-exchange amplitude could be comparable in size to the (Ge p ) 2 term in the unpolarized cross section 
at large Q 2 . In contrast, it was found that the two-photon exchange effects do not impact the polarization- 
transfer extraction of Ge p /G p m in an equally significant way. Thus a missing and un-factorisable part of the 
two-photon exchange amplitude at the level of a few percent may well explain the discrepancy between the 
two methods. 

Realistic calculations of elastic electron-nucleon scattering beyond the Born approximation are required in 
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order to demonstrate in a quantitative way that 27 exchange effects are indeed able to resolve this discrepancy. 

Recently, several model calculations of the 27 exchange amplitude have been done. In [Blu03|, a cal- 
culation of the 27 exchange when the hadronic intermediate state is a nucleon was performed. It found that 
the 27 exchange correction with intermediate nucleon can partially resolve the discrepancy between the two 
experimental techniques. However, subsequently it was found in [Kon05] that the effect is partly canceled 
when including the next hadronic intermediate state, the A(1232) resonance. The 27 exchange contribution to 
elastic eN scattering has also been estimated at large momentum transfer [Che04, Afa05a], through the scat- 
tering off a parton in a proton by relating the process on the nucleon to the generalized parton distributions. 
This approach effectively sums all possible intermediate states corresponding to excitations of the nucleon . 
Applying the two-photon exchange corrections to the unpolarized data (see dotted curve in Fig. [28]), yields 
a much flatter slope for the Rosenbluth plot, hence a much smaller value of Ge- The two-photon exchange 
corrections to the Rosenbluth process can therefore substantially reconcile the two ways of measuring Ge/ Gm 
(compare dotted with thin solid curves in Fig.l28l). 

To push the precision frontier further in electron scattering, one needs a good understanding, of 27 ex- 
change mechanisms, and of how they may or may not affect different observables. This justifies a systematic 
study of such 27 exchange effects, both theoretically and experimentally. Experimentally, the real part of the 
27 exchange amplitude can be accessed through the difference between elastic electron and positron scattering 
off a nucleon. Such experiments are planned in the near future. 

To conclude, on the one hand the discussion above makes it clear that the radiative corrections, including 
two hard photon exchange, for the cross section data are not complete at this point in time. Therefore, the FFs 
G\ and G 2 M obtained using the Rosenbluth method above Q 2 of 2 GeV 2 are not correct. On the other hand, 
all the authors cited above agree that radiative corrections change the longitudinal and transverse polarization 
components, P t and P e , in e + p — > e + p, similarly, with the ratio Pt/Pe affected only at the level of a 
few percent. The radiative corrections specifically calculated for the JLab polarization data by Afanasev et 
al. HAfaOll found that the corrections are ~1%, whereas the hard two photon exchange effects are at the few 
percent level HChe04[ Afa05aJ. Hence the polarization transfer method gives correct values for the FFs. 
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4 Theoretical interpretation of nucleon electromagnetic form factors 



In this section we give an overview of the theoretical understanding of the nucleon e.m. FFs. These FFs en- 
code the information on the structure of a strongly interacting many-body system of quarks and gluons, such 
as the nucleon. This field has a long history and many theoretical attempts have been made to understand 
the nucleon FFs. This reflects the fact that a direct calculation of nucleon FFs from the underlying theory, 
Quantum Chromodynamics (QCD), is complicated as it requires, in the few GeV momentum transfer region, 
non-perturbative methods. Hence, in practice it involves approximations which often have a limited range of 
applicability. Despite their approximations and limitations, some of these non-perturbative methods do reveal 
some insight in the nucleon structure. 

The earliest models to explain the global features of the nucleon FFs, such as its approximate dipole behav- 
ior, were vector meson dominance (VMD) models which are discussed in Sect. 14.11 In this picture the photon 
couples to the nucleon through the exchange of vector mesons. Such VMD models are a special case of more 
general dispersion relation fits, which allow to relate time-like and space-like FFs, and which are discussed 
subsequently. 

To understand the structure of the nucleon in terms of quark and gluon degrees of freedom, constituent 
quark models have a long history. We discuss the intricacies in describing a bound system of relativistic con- 
stituent quarks and review the resulting predictions for FFs in Sect. 14.21 Despite some of their successes, 
models based on quarks alone do suffer from the evident shortcoming that they do not satisfy the global chiral 
symmetry of QCD when rotating left and right handed light quarks in flavor space. This chiral symmetry 
is broken spontaneously in nature, and the resulting Goldstone bosons are pions. Since they are the lightest 
hadrons, they dominate the low momentum transfer behavior of form factors, and manifest themselves in a 
pion cloud surrounding the nucleon. Such pion cloud models will also be discussed in Sect. 14.21 

In Sect. 14.31 we discuss the spatial information which can be obtained from the nucleon FFs, and discuss 
both radial densities and the issue of shape of the nucleon. 

Sect. 14.41 describes the chiral effective field theory of QCD and their predictions for nucleon FFs at low 
momentum transfers, where such perturbative expansions are applicable. 

In Sect. 14.51 we shall discuss the lattice QCD simulations, which have the potential to calculate nucleon 
FFs from first principles. This is a rapidly developing field and important progress has been made in the recent 
past. Nevertheless, the lattice calculations are at present still severely limited by available computing power 
and in practice are performed for quark masses sizably larger than their values in nature. We will discuss the 
issues in such calculations and compare recent results. It will also be discussed how the chiral effective field 
theory can be useful in extrapolating present lattice QCD calculations to the physical pion mass. 

In Sect. 14.61 we discuss the quark structure of the nucleon and discuss generalized parton distributions 
(GPDs) of the nucleon. These GPDs are being accessed in hard exclusive reactions, which allow to remove 
in a controlled way a quark from the initial nucleon and implanting instead another quark in the final nucleon. 
The resulting GPDs can be interpreted as quark correlation functions and have the property that their first mo- 
ments exactly coincide with the nucleon FFs. We discuss the information which has been obtained on GPDs 
from fits of their first moments to the precise FF data set. 

Finally, in Sect. 14.71 we discuss the nucleon FFs in the framework of perturbative QCD. These consider- 
ations are only valid at very small distances, where quarks nearly do not interact. In this limit, the nucleon 
FFs correspond to a hard photon which hits a valence quark in the nucleon, which then shares the momentum 
with the other (near collinear) valence quarks through gluon exchange. We discuss the predictions made in 
this limit and confront them with the experimental status for Dirac and Pauli FFs at large momentum transfers. 
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4.1 Dispersion theory 



4.1.1 Vector Meson Dominance (VMD) 

The starting point in understanding the interaction of a vector probe such as the photon with a hadronic system 
is provided by the observation that the lowest lying hadrons with vector quantum numbers are the vector 
mesons p(770), c<j(782) and 0(1020). In the process e + e _ — > hadrons, these vector mesons show up as 
prominent resonances at the corresponding values of the e + e~ squared cm. energy q 2 > 0. One therefore 
expects that in the elastic electron scattering process on the nucleon, eiV — > eiV, the nucleon electromagnetic 
FFs at low space-like momentum transfers, q 2 < 0, will be dominated by these lowest lying singularities from 
the time-like region. A large class of models for Fx and F 2 are based on this vector meson dominance (VMD) 
hypothesis, as depicted in Fig. [291 



p, u, 4' 



Figure 29: The vector meson dominance picture for the coupling of the photon (with four-momentum q) to a nucleon. 

Within such VMD models, the approximate dipole behavior of the nucleon e.m. FFs, see Eq. (fl4l) . can be 
understood as being due to the contribution of two nearby vector meson poles which have opposite residua. 
Assume that one considers two vector meson pole contributions in Fig. [29] (with masses myi and my 2 and 
residua of equal magnitude and opposite sign a and —a respectively), one obtains : 



F 2 ( Q 2 ) - + = ^y""vi-""v2j (35) 
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An early VMD fit was performed by Iachello et al. [|Iac73ll and predicted a linear decrease of the proton 
Ge p /Gm p ratio, which is in basic agreement with the result from the polarization transfer experiments. Such 
VMD models have been extended by Gari and Kriimpelmann [Gar85] to include the perturbative QCD (pQCD) 
scaling relations [ |Lep80[ , which state that (see Sect. [477]) F x ~ 1 /Q 4 , and F 2 /F 1 ~ 1 /Q 2 . 



In more recent years, extended VMD fits which provide a relatively good parameterization of all nu- 
cleon e.m. FFs have been obtained. An example is Lomon's fit HLomOlll , using p(770), u(782), 0(1020), 
and p'(1450) mesons and containing 11 parameters. Another such recent parameterization by Bijker and 



Iachello [Bij04| including p(770), u(782), and 0(1020) mesons only achieves a good fit by adding a phe- 
nomenological contribution attributed to a quarklike intrinsic qqq structure (of rms radius ~ 0.34 fm) besides 
the vector-meson exchange terms. The pQCD scaling relations are built into this fit which has 6 free parame- 
ters which are fit to the data. In contrast to the early fit of Ref. HIac73H . the new fit of Ref. [Bij04| gives a very 
good description of the neutron data, albeit at the expense of a slightly worse fit for the proton data. 

It will be interesting to check the resulting VMD fits for the neutron FFs to larger Q 2 . In this regard, an 
interesting "prediction" can be drawn when the FFs F 2 and Fx obtained directly from double polarization ex- 
periments are shown in the same graph for the proton and the neutron, as in Fig. [30] It is remarkable that both 
Fx and F 2 tend toward the same value for proton and neutron, and may meet at a Q 2 value which will soon be 
accessible for the neutron. This conclusion is influenced by the VMD fits shown in the same figure, and rests 
on their extrapolation for the neutron to larger Q 2 . Note that the VMD fits shown include all data for p and n, 



but selects the recoil polarization over the Rosenbluth results for Q 2 larger than 1 GeV 2 
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Figure 30: The FFs F\ and F2 for the proton and the neutron obtained from double polarization experiments only. The 
values of F2 and F\ were obtained from the experimental FF ratios using fitted values to the data base for Gmp and Gmu- 
For the proton data, the fit from [Kel04] was used; for the neutron data in [Mad03] the fit from [Kel04], and in [War04] 
the fit from HKub021 were used. The curves are the VMD fits of Lomon [LomOl ] and of Bijker and Iachello |Bij04 1. 



4.1.2 Dispersion analyses 

Despite the relatively good fits obtained by the VMD models, it was already pointed out in 1959 by Frazer 
and Fulco HFra59l that such an approach is at odds with general constraints from unitarity. Assuming an 
unsubtracted dispersion relation (DR), the nucleon e.m. FFs F(q 2 ), where F generically stands for any of the 
four FFs, can be obtained as : 

, 9s 1 [°° , , 9 ImF(g' 2 ) 

The dispersion analyses are performed separately for nucleon isoscalar and isovector FFs. In the vector- 
isovector spectral function ImF(g' 2 ) one notices a large non-resonant contribution starting from t = 4m 2 
and extending under the p-peak. Such a non-resonant contribution arises due to the two-pion continuum. For 
the isoscalar spectral function, the integral starts at t = 9m 2 , corresponding to 3ir intermediate states. The 
two-pion continuum contribution was estimated by Hohler and collaborators [Hoh76J by using pion time-like 
FF data and nn — > NN amplitudes which were determined by extrapolating nN partial waves to the time-like 
region [Hoh75l. 

Holder's analysis has been updated by Mergell, Meissner, and Drechsel [Mer96| in the mid-nineties and 
extended to include the nucleon time-like FF data [Ham96 |. The inclusion of recent neutron FF data in such 
dispersion relation analysis has been performed in Ref. [Ham04]. The resulting analysis describes the nucleon 
isovector FFs through the 2n continuum (including the p(770)), and three additional vector isovector meson 
poles : p'(1050), p"(1465), p'"(1700). The isoscalar FFs are described by four vector isoscalar meson poles : 
w(770), 0(1020), S"(1650) and S"'(1680). In this approach, the masses of the mesons p', p", p'" , S\ S" and 
the 14 residua (one for both the vector (F\) and tensor (F 2 ) channels for each meson) are fitted and the pQCD 
scaling behavior is parameterized through three additional parameters. Note that for the isovector channel, the 
fitted masses for p" and p'" correspond to physical particles listed by the Particle Data Group (PDG), whereas 
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enforcing the correct normalization of all FFs, the experimental value for the neutron charge radius, as well 
as the pQCD scaling behavior, requires the inclusion of an unphysical p' meson with mass 1050 MeV. The 
analysis of Hammer and Meissner [Ham04] also finds that the residua for both isovector FFs FY and F 2 y of p" 
and p'" are relatively close in magnitude and opposite in sign, required by the approximate dipole behavior of 
the isovector FFs. For the isoscalar FFs Ff and F 2 5 , the fit also drives the residua of the nearby poles S' and 
S" to values very close in magnitude and of opposite signs, required by the approximate dipole behavior of the 
isoscalar FFs. Using such an analysis, a good description of most FF data with the exception of the Ge p /Gm p 
polarization data at Q 2 > 3 GeV 2 , was obtained in [Ham04l. 
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Figure 31: Dispersion relation (15 parameter) fit for the four nucleon (space-like) e.m. FFs compared with the world 
data (circles) including the JLab/CLAS data for Gun (triangles) [Bro05]. The dashed curves indicate the la deviation 
from the fit, given by the solid curves. Figure from Ref. HBel07ll . 



The dispersion relation analysis of nucleon e.m. FFs has been further improved by Belushkin et al. HBel071l . 
In addition to the 2n continuum present in the isovector spectral functions of the previous DR analyses, also the 
pn and KK continua were included as independent input in the isoscalar spectral functions. In Ref. [B el07H . 
the 2tt continuum was reevaluated using the latest experimental data for the pion FFs in the time-like re- 
gion [Bel06J . The KK continuum was obtained from an analytic continuation of KN scattering data [Ham99 1 . 
Following the work of Ref. [Mei97], the pit continuum was approximated in the DR analysis by an effective 
pole term for a fictitious to' meson with mass : m w / = 1.12 GeV. This approximate pn continuum is found 
to yield an important negative contribution to Ff . The remaining contributions to the spectral functions are 
parameterized by vector meson poles from a fit to the FF data. The parameters in the fit were constrained 
to yield the correct normalization of the FFs at zero momentum transfers. The asymptotic constraints from 
pQCD were included in two different forms : either as a superconvergence relation or by adding an explicit 
continuum term with the imposed pQCD behavior. A simultaneous fit to the world data for all four FFs in 
both the space-like and time-like regions was performed. Fig. |3T] shows this fit for the nucleon space-like FFs 
where for Ge p /Gm p at larger Q 2 the JLab/Hall A polarization data HJonOOl |Gay02[ |Pun 05 1 have been used, 
and where for Gun the preliminary JLab/CLAS data [Bro05] have been included. In this fit, the pQCD limit 
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was imposed through an explicit continuum term and the minimum number of poles in addition to the irn, pn 
and KK continua were chosen to fit the data. In addition to the u;(782), the fit yields two more isoscalar poles 
(m s i ~ 1.05 GeV and m s2 — 1.4 GeV), and three additional isovector poles (m v i ~ 1.0 GeV, m v2 — 1.6 GeV, 
and m V 3 ~ 1.8 GeV ). The resulting 15 parameter fit shown in Fig.[3T]has a total x 2 /d.o.f. value of 2.2. 

It will be interesting to confront the most recent and sophisticated DR fit of [Bel07J with upcoming data for 
G E p/G M p out to 8.5 GeV 2 HBra04H . In all discussed VMD and DR fits starting with Gari-Krumpelmann HGar85l 
the asymptotic pQCD limit F 2 /Fi ~ l/Q 2 was built in, although the data do not support this limit at avail- 
able momentum transfers, see Fig.|43]in Sect. 14.71 Besides, there is various theoretical work indicating that 
the pQCD prediction, in particular for F 2p , might only set in at significantly larger values of Q 2 , of the order 
of several tens of GeV 2 . It might therefore be worthwhile to investigate how the DR analysis changes by 
removing the bias F 2 jF\ ~ l/Q 2 from the analysis, when fitting data in the range up to Q 2 ~ 10 GeV 2 . 

4.2 Quark models versus pion-cloud models 
4.2.1 Constituent quark models 

In our quest to understand the structure of the nucleon in terms of the quark and gluon degrees of freedom 
which appear in the QCD Lagrangian, constituent quark models (CQMs) have a long history, which predates 
the establishment of the theory of strong interactions, QCD. In a CQM, the nucleon appears as the ground 
state of a quantum-mechanical three-quark system in a confining potential. In such a picture, the ground state 
baryons (composed of the light up (it), down (d) and strange (s) quark flavors) are described by SU (6) spin- 
flavor wave functions (WFs), supplemented by an antisymmetric color WF. 



In the Isgur-Karl model [Isg78 1, the constituent quarks move in a harmonic oscillator type confining poten- 
tial. For the ground state baryons, the three constituent quarks are in the Is oscillator ground state, correspond- 
ing to the [56]-plet of SU (6). In the Isgur-Karl model, the long-range confining potential is supplemented by 
an interquark force corresponding to one-gluon exchange. The one-gluon exchange leads to a color hyper- 
fine interaction between quarks, which breaks the SU(6) symmetry and leads to a mass splitting between 
iV(939) and A (1232), often referred to as the hyperfine splitting. It was found that it also predicts well the 
mass splittings between octet and decuplet baryons [DeR75|. Furthermore, the color hyperfine interaction 
leads to a tensor force which produces a small /J-state (L = 2) admixture in the N (as well as A) ground 
states [Kon80, Isg82[, corresponding to a TJ-state probability in the N ground state around 0.2 %. Even 



though such TJ-wave probability is small, it leads to a non-spherical charge distribution. For a static charge 
distribution, a measure of the non- sphericity (or deformation) is given by its quadrupole moment. Since the 
nucleon has spin 1/2, an intrinsic quadrupole moment of the nucleon cannot be directly measured because 
angular momentum conservation forbids a non-zero matrix element of a (L = 2) quadrupole operator between 
spin 1/2 states. However this quadrupole deformation may reveal itself in an electromagnetically induced tran- 
sition from the spin 1/2 N to the spin 3/2 A state. In this way, the tensor force between quarks gives rise to 
non-zero values for the electric quadrupole (E2) and Coulomb quadrupole (C2) transitions^. 

The non-relativistic CQM, despite its simplicity, is quite successful in predicting the spectrum of low-lying 
baryons, and gives a relatively good description of static properties such as the octet baryon magnetic mo- 
ments. To calculate the FFs of a system of constituents with masses small compared with the confinement 
mass scale necessitates however a relativistic treatment even for low momentum transfers. For momentum 



3 The relation between the tensor force, Z?-wave admixture, and the electromagnetic N — > A transition was already pointed out 
in the early paper of Glashow |Gla79|. An up-to-date discussion of this field can be found in the review of |Pas07|. 
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transfers several times the nucleon mass squared a relativistic description becomes crucial. 

In contrast to the calculation of the spectrum, which uses eigenfunctions of a Poincare invariant mass op- 
erator, a calculation of the nucleon electromagnetic FFs requires the relation between the rest frame spin and 
momenta (in the three-quark WF) and those in the moving frame. This requires an extension of eigenfunc- 
tions of the spin and mass operators, so as to transform consistently under the unitary representations of the 
Poincare group. The way to implement relativity into a Hamiltonian formalism (describing e.g. a system of 
three interacting constituent quarks) has been laid out by Dirac [Dir49]. There are three forms of dynamics 
(so called instant, point, and light-front forms) which differ in the choice of the kinematical subgroup of the 
Poincare group. This is the subgroup of the Poincare group whose commutator relations are not affected by 
the interactions between the constituents. The three (unitarily equivalent) forms therefore differ by which of 
the ten generators of the Poincare group (four space-time translations, three spatial rotations, and three boosts) 
are kinematical (i.e. interaction free), and which are dynamical, i.e. depend on the interactions and necessarily 
have to be approximated in a practical calculation. 

In the instant form, the dynamical generators are the time component of the four-momentum and the three 
boost operators. Rotations do not contain interactions, which makes it easy to construct states of definite an- 
gular momentum in this form. 

In the point form, both boosts and rotations are kinematical. The point-form therefore has the important 
technical advantage that the angular momenta and Lorentz boosts are the same as in the free case. However all 
four components of the four- vector operator are dynamical in this form. 

In the light-front form, seven of the generators of the Poincare group are kinematical (this corresponds to 
the symmetry group of a null plane), which is the maximum number possible. The remaining three dynamical 
generators which contain the interactions are one component of the four-momentum operator (the so-called 
light-cone Hamiltonian) and 2 transverse rotations. Light-front (as well as point form) calculations for rela- 
tivistic CQMs are convenient as they allow to boost quark WFs independently of the details of the interaction. 
The drawback of the light- front calculations however is that because two generators of rotations are dynamical, 
the construction of states with good total angular momentum becomes interaction dependent. 

Any practical calculation in one of the three forms approximates the current operator. The common (so- 
called impulse) approximation is that the photon interacts with a single quark in the nucleon. 

The light-front form calculation of nucleon FFs has been pioneered by Berestetsky and Terentev [B er76H . 
and more recently developed by Chung and Coester HChu91H . In practice one starts from a rest frame nucleon 
WF for the three-quark state which ideally is fitted to the baryon spectrum. The nucleon WF in the light-front 
form (so-called light-front WF) is obtained by a Melosh rotation [Mel74| of each of the quark spinors, con- 
necting the instant and light-front forms. When performing the front form calculation in a (Drell-Yan) frame 
where the photon light-cone momentum component q + = 0, the space-like virtual photon only connects 
Fock components in the nucleon light-front WFs with the same number of constituents, i.e. matrix elements 
between qqq and qqqqq states which would be present in an instant form calculation are zero in the light-front 
calculation. This property allows for a consistent calculation within the light-front formalism when truncating 
the Fock space to only the three-quark state. 

In HChu91H a Gaussian WF in the quark internal (transverse) momentum variables was used. Although 
this model yields a surprisingly good agreement for the observed Ge p /Gm p ratio, see Fig. [32l it yields nu- 
cleon FFs which drop too fast at larger Q 2 values when using constituent quark masses around 330 MeV. 
Schlumpf [Schl93| allowed for high momentum components in the nucleon light-front WF by adopting a 



Defining light-cone components as x ± = (x° ± x 3 )/\/2 and defining the null-plane by x + = 0. 
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power law dependence in the quadratic quark internal momentum variables. The two parameters in Schlumpf 's 
WF were fitted to magnetic moments and semi-leptonic decays of the baryon octet. The resulting e.m. FF cal- 
culations reproduce reasonably well the power behavior of the FF at larger Q 2 . The WF of Schlumpf was 
also used by Frank, Jennings, and Miller [Fra96, Mil02aJ. They found that such a light-front WF leads to a 
violation of hadron helicity conservation resulting in a F 2p / F\ p ratio which drops less fast than 1/Q 2 [Mil02a], 
in agreement with the G Ep /G Mp polarization data. 
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Figure 32: Comparison of relativistic CQM calculations with the data for h p Ge p /Gm p - Dotted curve : front form 
calculation of Chung and Coester [Chu91] with point-like constituent quarks; thick solid curve : front form calculation 
of Frank et al. |Fra96]; dot-dashed curve : front form calculation of Cardarelli et al. MCar95[ ICarOOl with point-like 
constituent quarks; dashed curve : point form calculation of Boffi et al. [BofOl ] in the Goldstone boson exchange model 
with point-like constituent quarks; thin solid curve : covariant spectator model of Gross and Agbakpe [Gro06|. The data 
are from BPun05l (solid circles) and [Gay02 ] (empty squares). 



The WFs in the calculations described above were however not constructed from a detailed fit to the baryon 
spectrum. Cardarelli et al. subsequently performed a more "microscopic" light-front calculation [Car95, 
ICarOOB where the light-front WF was obtained from a rest frame WF which provided a fit to the spectrum. The 



rest frame WF was taken from the relativized Capstick-Isgur model [Cap86|. Using this WF, the constituent 
quark momentum distribution in the nucleon was found to yield an important content of high-momentum com- 
ponents, which are generated by the short-range part of the quark-quark interaction, which is due to one-gluon 
exchange in the Capstick-Isgur model. These components are completely absent if one only considers the 
linear confinement potential in the model. 

In a CQM calculation, the effect of other degrees of freedom beyond three quarks are buried within the 
constituent quarks, which are considered as quasi-particles. In the absence of a microscopic calculation, such 
effects are parameterized in terms of constituent quark FFs. In HPet031l . it was shown that the data for the 
proton unpolarized forward structure function at low momentum transfers exhibits a new scaling property and 
can be interpreted as quasi-elastic scattering off extended constituent quarks inside the proton described by a 
constituent quark FF. The resulting constituent size is around 0.2 - 0.3 fm. Using such effective constituent 
quark FF in the light-front form calculation of HCarOOL allows a good description of the individual nucleon 
FFs, see [PacOOJ . Note however that the experimental G e p / G m p ratio can basically be reproduced using point- 
like constituent quarks, see Fig. [32] The suppression of the Ge p /Gm p ratio with respect to the dipole-fit as 
predicted in the light-front form CQM calculation is attributed to relativistic effects generated by the Melosh 
rotations of the constituent quark spins. These Melosh rotations introduce kinematical SU (6) breaking effects 
in addition to the dynamical SU (6) breaking due to the (hyperfine) one-gluon exchange potential. 
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A comparable amount of high-momentum components in the nucleon WF was obtained in the Goldstone - 
boson-exchange (GBE) quark model [Glo98a, Glo98b j. This model relies on constituent quarks and Goldstone 
bosons, which arise as effective degrees of freedom of low-energy QCD from the spontaneous breaking of the 
chiral symmetry. The resulting CQM assumes a linear confinement potential supplemented by a quark-quark 
interaction based on the exchange of pseudoscalar Goldstone bosons, which is the source of the hyperfine 
interaction. It was shown in [Glo98a, Glo98b] that the GBE CQM yields a unified description of light- and 
strange-baryon spectra. The GBE CQM was used in [Wag011 [Bof01|| to calculate the nucleon e.m. FFs in the 
point-form. The neutron charge radius is well described in this model and is driven by the mixed-symmetry 
component in the neutron WF. In contrast to the light-front calculation [CarOO, PacOOJ, it was found that when 
performing a point-form calculation of the nucleon e.m. FFs at larger Q 2 within the impulse approximation, 
i.e. considering only single-quark currents, a surprisingly good overall description of the nucleon e.m. FFs 
can be obtained, using point-like constituent quarks only. When looking at details of Refs. [WagOl, BofOl], 
the agreement is worse though for Gm p which is underpredicted at larger Q 2 , and the ratio of Ge p /Gm p is 
overpredicted at larger Q 2 , see Fig. |32] Similar findings have also been obtained in the point-form calculation 



of [ Wag05 1 for the OGE CQM. The overall success of the point-form result using point-like constituent quarks 
was attributed in [ WagO fj IBofO 1 [ |Wag05[ to the major role played by relativity. Such a finding is remarkable 
in view of the expected finite size of the constituent quarks, as discussed above. 

An explanation for the above finding for the nucleon e.m. FFs in the point form, using the single-quark 
current approximation, has been suggested by Coester and Riska HCoe03L When the spatial extent of the 
three-quark WF is scaled (unitarily) to zero, both instant- and front-forms yield FFs independent of the mo- 
mentum transfer. Therefore, to reproduce the experimental fall-off of the nucleon e.m. FFs at large momentum 
transfers requires the introduction of constituent quark FFs. In contrast, when the WF in point form is scaled 
unitarily to zero (so-called point limit), a non-trivial scaling limit is obtained for the FFs, depending on the 
shape of the WF. At high values of momentum transfer, the scaled FFs decrease with an inverse power of the 
momentum transfer. The power is determined by the current operator and is independent of the WF. An explicit 
comparative calculation of the baryon e.m. FFs between the three different forms was performed in [Jul04J 
using a simple algebraic form for the three-quark WF, depending on two parameters. It was verified that a 
qualitative description of the nucleon FF data demands a spatially extended WF in the instant- and front-form 
descriptions, in contrast to the point-form description which demands a much more compact WF. 

A manifestly covariant CQM calculation within the Bethe-Salpeter formalism and using an instanton- 
induced interaction between quarks has been performed by Merten et al. [Mer02J. Although this model repro- 
duces the baryon spectrum, it can only qualitatively account for the Q 2 dependence of the nucleon e.m. FFs. 

Another covariant CQM calculation was performed by Gross and Agbakpe [Gro06|, using a covariant spec- 
tator model. Assuming a simple pure S-wave form for the nucleon three-quark wave function, evaluating the 
current matrix element in a relativistic impulse approximation, and assuming constituent quark FFs including 
a phenomenological term which parameterizes the pion cloud, an eleven parameter description of the nucleon 
FF data was obtained, see Fig. [32] 

As a next step for CQMs, it would clearly be very worthwhile to investigate the approximations in the cur- 
rent operator within each form. The quality of the commonly made impulse approximation may differ between 
the different forms. Within the context of a toy model calculation in Refs. [De s04[ IDes06L it has e.g. been 
shown that the neglect of two-body currents in the point form does affect the FFs in a more drastic way than 
their neglect in the instant or light-front forms. 

The importance of two-body currents has also been shown in the work of De S anctis et al. [ DeSOO | . In that 
work, a calculation within the hypercentral CQM was performed of the (two-body) quark pair contribution 
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to the e.m. current resulting from the one-gluon exchange interaction between the quarks. This pair current 
contribution was found to lead to a sizeable reduction of Ge p compared with Gm p - 



4.2.2 Pion cloud models 

Despite their relative success in describing the spectrum and structure of low-lying baryons, models based on 
constituent quarks alone suffer from evident shortcomings as they do not satisfy all symmetry properties of the 
QCD Lagrangian. In nature, the up and down (current) quarks are nearly massless. In the exact massless limit, 
the QCD Lagrangian is invariant under SU(2) L x SU(2) R rotations of left (L) and right (R) handed quarks 
in flavor space. This chiral symmetry is spontaneously broken in nature leading to the appearance of massless 
Goldstone modes. For two flavors, there are three Goldstone bosons — pions, which acquire a mass due to the 
explicit breaking of chiral symmetry by the current quark masses. 

Since pions are the lightest hadrons, they dominate the long-distance behavior of hadron WFs and yield 
characteristic signatures in the low-momentum transfer behavior of hadronic FFs. Therefore, a natural way to 
qualitatively improve on the above-mentioned CQMs is to include the pionic degrees of freedom [Man84j. 

An early quark model with chiral symmetry is the chiral (or, cloudy) bag model. This model improves the 
early MIT bag model by introducing an elementary, perturbative pion which couples to quarks in the bag in 
such a way that chiral symmetry is restored [Tho82]. Within the cloudy bag model, Lu et al. flLu98H performed 
a calculation of the nucleon e.m. FFs improving upon previous calculations by applying a correction for the 
center-of-mass motion of the bag. This calculation also implemented Lorentz covariance in an approximate 
way by using a prescription for the Lorentz contraction of the internal structure of the nucleon. Using a bag 
radius R ~ 1 fm, this model provides a good description of the nucleon e.m. FFs in the range Q 2 < 1 GeV 2 . 

To extend such a calculation to larger Q 2 , Miller performed a light-front cloudy bag model calcula- 
tion [Mil02bj. Starting from a model in terms of constituent quarks [Mil02a], described by the light-front 
WF of Schlumpf, the effects of the pion cloud were calculated through one-loop diagrams, including relativis- 
tic nNN vertex FFs. The model gives a relatively good gobal account of the data both at low Q 2 and larger 
Q 2 , though tends to show too much structure around the dipole form for the magnetic FFs at low Q 2 . 

The cloudy bag model is one chiral quark model which treats the effect of pions perturbatively. Other 
quark models which calculated nucleon e.m. FFs using perturbative pions can be found e.g. in the early works 
of HOse84[ |Jen92L as well as in the already discussed works of [Glo98a, Glo98b|. Recently, the above chiral 
quark models where pions are included perturbatively have been improved in [Fae06a]. This work extends a 



previous work of [LyuOl | by dynamically dressing bare constituent quarks by mesons to fourth order within a 
manifestly Lorentz covariant formalism. Once the nucleon and A hyperon magnetic moments are fitted, other 
e.m. properties, such as the nucleon e.m. FFs at low momentum transfers, follow as a prediction. It was found 
in HFae06all that the meson cloud is able to nicely describe the FF data in the momentum transfer region up to 
about 0.5 GeV 2 . To extend the calculations to larger Q 2 , a phenomenological approach has been adopted in 
[Fae06a] by introducing bare constituent quark FFs which were parameterized in terms of 10 parameters. Such 
parameterization makes it plausible to simultaneously explain the underlying dipole structure in the nucleon 
e.m. FF as well as the meson cloud contribution at low Q 2 which results from the underlying chiral dynamics. 
In a later paper [Fae06b|, a model calculation for the bare constituent quark FFs has been performed and ap- 
plied to the e.m. properties of the N — ► A transition. 

When pion effects dominate nucleon structure, their effects have to be treated non-perturbatively. A non- 
perturbative approach which has both quark and pion degrees of freedom and interpolates between a CQM and 
the Skyrme model (where the nucleon appears as a soliton solution of an effective nonlinear pion field theory) 
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is the chiral quark soliton model (xQSM). As for the Skyrme model, the xQSM is based on a l/N c expansion 
(with N c the number of colors in QCD). Its effective chiral action has been derived from the instanton model 
of the QCD vacuum HDia86il . which provides a natural mechanism of chiral symmetry breaking and enables 
one to generate dynamically the constituent quark mass. Although in reality the number of colors N c is equal 
to three, the extreme limit of large N c is known to yield useful insights. At large N c the nucleon is heavy and 
can be viewed as N c "valence" quarks bound by a self-consistent pion field (the "soliton") HDia 88|. A suc- 
cessful description of static properties of baryons, such as mass splittings, axial constants, magnetic moments, 
FFs, has been achieved (typically at the 30 % level or better, see [Chr96] for a review of early results). After 
reproducing masses and decay constants in the mesonic sector, the only free parameter left to be fixed in the 
baryonic sector is the constituent quark mass. When taking rotational (l/N c ) corrections into account, this 
model achieved a qualitative good description of the nucleon e.m. FFs in the range Q 2 < 1 GeV 2 , using a 
constituent quark mass around 420 MeV HChr95L The chiral soliton model naturally accounts for the decrease 
of the Gep/Gmp ratio with increasing Q 2 . This can be understood from the hedgehog structure in soliton 
models which couples spatial rotations with isorotations. In the soliton rest frame, the iso vector electric FF 
G\ therefore measures the rotational inertia density p v (r), in contrast to the isoscalar electric FF G S E which 
measures the isoscalar baryon density p s {r). For a rigid rotor, the inertia density is obtained from the baryon 
density as p v = r 2 /r 2 B p s , with r B a free parameter characterizing the spatial extent. Assuming a Gaussian 
density for p s (r), this yields [Hol96| : 

G Mp (Q 2 ) - 1 18 QrB - (37) 

With the choice r 2 B « (0.3 fm) 2 , one can obtain an excellent fit of the polarization data for Ge p /Gm p - Although 
in the chiral soliton model calculation the baryon density is not exactly Gaussian, and the rigid rotor calculation 
does not hold exactly, these relations can be considered as approximate relations [Hol96]. 

Holzwarth [Hol96J extended the chiral soliton model by including the p (to) meson propagators for the 
isovector (isoscalar) channels respectively. Furthermore, to extend the range in Q 2 of the predictions, he uses 
a relativistic prescription to boost the soliton rest frame densities to the Breit frame. Such prescription is also 
used to extract radial charge and magnetization rest frame densities from experimental FFs, as will be discussed 
in Sect. 14.31 Using 4 fit parameters (one effective boost mass and three free parameters to fix the couplings of 
p and to mesons), the model was found to provide a good account of the detailed structure of the nucleon e.m. 
FFs in the low Q 2 region. In particular, for Ge p /Gm p it predicts a decreasing ratio in good agreement with 
the data. At larger Q 2 , the boost prescription gives a reasonably good account of the data (except for G Mn ) 
and predicts a zero in Ge p around 10 GeV 2 . Due to the uncertainty introduced from the particular choice for 
the boost prescription, the high Q 2 behavior (for Q 2 larger than about 4M 2 ) of the e.m. FFs is however not a 
profound prediction of the low-energy effective model. 



4.3 Radial distributions and shape of the nucleon 

As discussed in Sect. 12.1.31 in the Breit frame the nucleon charge operator depends only on the electric FF 
Ge, whereas the e.m. three-current operator depends only on the magnetic FF Gm- This suggest to interpret 
the Fourier transforms of Ge (Gm) as the nucleon charge (magnetization) densities. This identification is only 
appropriate for a non-relativistic (static) system however, as in general there is a variation of the Breit frame 
with Q 2 . For the nucleon, where FF data have been obtained for Q 2 values much larger than M 2 , one needs 
to take the effect of relativity into account. Recently Kelly [Kel02] has used a relativistic prescription to relate 
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the Sachs FFs to the nucleon charge and magnetization densities, accounting for the Lorentz contraction of the 
densities in the Breit frame relative to the rest frame. 

One starts from the spherical nucleon charge p c h(r) and magnetization p m (r) densities in the nucleon rest 
frame. These densities are normalized so as to yield the total charge for p c h, or one for p m (the magnetic 
moment is taken out of the density) as : 

POO 

drr 2 p ch (r) = Z, / drr 2 p m (r) = 1, (38) 

o Jo 

where Z — 0, 1 is the nucleon charge. From these intrinsic (rest frame) densities, one can construct intrinsic 
FFs p{k) which are related through a Fourier-Bessel transform as : 



p(k)= drr 2 j (kr)p(r), p{r) = - / dk k 2 j Q (kr)p(k), 

Jo 71 Jo 



(39) 



with k = |g| is the wave vector in the nucleon rest frame. For a non-relativistic system, the intrinsic FFs are 
obtained from the Sachs FFs using k — > Q as : p c h(k) — > G E (Q 2 ), and pNpm(k) — > Gm(Q 2 )- 

To properly relate the intrinsic FFs evaluated in the rest frame to the Breit frame, where the nucleon moves 
with velocity v = a/t/(1 + t) relative to the rest frame, involves a Lorentz boost with : 7 2 = (1 — v 2 )^ 1 = 1 + 
r. This Lorentz boost leads to a contraction of the nucleon densities as seen in the Breit frame. Consequently, 
in the Fourier-transforms, this amounts to replace in the intrinsic FF arguments : 

k 2 ^Q 2 /(l + r). (40) 

To relate intrinsic FFs p(k) with the Sachs FFs G(Q 2 ) is not unambiguous however because the boost operator 
for a composite system depends on the interactions among its constituents. There exist different prescription 
in the literature which can be written in the form : 

Pch(k) = 1 2nE G E (Q 2 ) = (l+T) nE G E (Q 2 ), (41) 
p N p m {k) = l 2nM G M (Q 2 ) = (l + T) n ^G M (Q 2 ), (42) 

where k and Q 2 are related as in Eq. (|40l ). For Q 2 — > oo, the boost maps G(Q 2 — > oo) to p(2M). One sees that 
there is a limiting wave vector k max = 2M determined by the nucleon Compton wavelength. In the rest frame, 
no information can be obtained on distance scales smaller than the Compton wavelength due to relativistic 
position fluctuations (known as the Zitterbewegung). To account for an asymptotic l/Q A FF behavior, Mitra 
and Kumari [MitWJ proposed the choice n E = um = 2. Kelly [Kel02J followed this choice when extracting 
the rest frame densities from the measured nucleon e.m. FFs. 

In his analysis, Kelly furthermore minimized the model dependence of the fitted densities by using an ex- 
pansion in a complete set of radial basis functions. For Q 2 > 1 GeV 2 the G Ep analysis used recoil polarization 
data from JLab [JonOOl |Gay02[ rather than the Rosenbluth separation data. Fig. [33] compares the fitted charge 



and magnetization densities for neutron and proton. The uncertainty bands include both statistical and incom- 
pleteness errors. The excess of negative charge near r ~ 0.8 — 1.0 fm is a characteristics of the 7r~-meson 
cloud in the neutron. The proton charge density is significantly broader than the magnetization density, a direct 
consequence of G Ep being softer than G Mp in Q 2 -space. 

To investigate the pion cloud as revealed through the neutron electric charge distribution further, Friedrich 
and Walcher [Fri03J have performed a phenomenological analysis of all four nucleon e.m. FFs. They per- 
formed a two-component fit of the four FFs starting from a smooth part (parameterized by a sum of two 
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Figure 33: Radial distributions of charge p c h and magnetization p m in the proton and neutron from [Kel02], obtained 
from the Fourier Bessel transforms Eq. d39l > using the relativistic transformation of Eqs. (14 11421 ) with tie = um = 2. 
Note that the neutron charge distribution has been multiplied by a factor of 6 to emphasize the similarity in shape of 
charge and magnetization densities. 

dipoles) and by adding on top of it a Gaussian "bump" structure. The choice of such a two-component form 
was triggered by the behavior of Geu at small Q 2 , and by the observation of the noticeable oscillations of the 
other three e.m. FFs around the dipole form. Their parameterization allows for 6 fit parameters for each FF, 
which provide an excellent fit to the FFs. When subtracting from the data two dipoles with suitably chosen 
parameters, the remaining part displays a bump structure as shown in Fig.[34l Friedrich and Walcher made the 
striking observation that all four FFs display such a bump structure around Q 2 w 0.25 GeV 2 . They intepret this 
structure as a signature of the pion cloud. Upon Fourier transforming, the corresponding Breit-frame densities, 
corresponding to the "bump" structure in the FFs, were found to extend as far out as 2 fm. It is interesting 
to compare this with the findings of the dispersion theory, in which the longest range part of the pion cloud 
contribution to the nucleon e.m. FFs is given by the 2n continuum. These 2n continuum contributions were 
found to be much more confined in coordinate space [Bel06J. In order to get a bump structure in Geu in the 
DR theory requires to introduce additional strength in the spectral functions below 1 GeV. New high precision 
data for Ge p /Gm p from the BLAST experiment at MIT-Bates HCra061 . shown in Fig. [22l confirms the dip 
structure around Q 2 fa 0.3 GeV 2 . It will also be interesting to compare upcoming data of BLAST for Geu m 



Miller [Mil03J has defined spin-dependent quark densities as matrix elements of density operators in proton 
states of definite spin-polarization. Within a constituent quark picture, the spin-dependent density operator for 
a quark in the proton to be found at position r and with spin-direction n is given by : 



where the sum runs over the three constituent quarks i with fractional charge ej/e. Relative to the spin-direction 
s of the proton, Miller then studied the distribution of quarks for different quark spin orientations n. The so 

5 For preliminary data from BLAST for Gb„, see Ref. |Zis05 1. 
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Figure 34: Phenomenological two-component fit of the nucleon e.m. FFs according to Friedrich and Walcher [Fri03|. 
The 6-parameter fit for each FF consists of a smooth part, described by a sum of two dipoles, and a Gaussian bump part. 
The latter is displayed in the figure. 



defined densities may become non-spherical as shown in Refs. [Mil 03llKvi06H . Averaging over quark spin n 
or over nucleon spin s yields a spherical distribution [Gro06J. 



4.4 Chiral perturbation theory 

At low momentum transfers Q 2 , the nucleon e.m. FFs can also be studied within chiral perturbation theory 
(xPT) expansions based on chiral Lagrangians with pion and nucleon fields. In xPT, the short-distance physics 
is parameterized in terms of low-energy-constants (LECs) which ideally can be determined by matching to 
QCD but in practice are fitted to experiment or are estimated using resonance saturation. In the calculation of 
the nucleon e.m. FFs, the LECs can be fitted to the nucleon charge radii and the anomalous magnetic moments. 
Once they are fixed, the Q 2 dependence of the FFs follows as a prediction. 

To calculate the nucleon e.m. FFs in %EFT involves a simultaneous expansion in soft scales : Q 2 and 
m n , which are understood to be small relative to the chiral symmetry breaking scale A x sb ~ 1 GeV. Several 
expansion schemes (also called power-counting schemes) have been developed in the literature. They all yield 
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the same non-analytic dependencies (e.g. terms proportional to m^, m 3 , m 2 him 2 , ...) but differ in analytic 
terms (e.g. terms proportional to m 2 , m\, ...). 

Because the first nucleon excitation, the A(1232) resonance, has an excitation energy of only about 
A = Ma — M ~ 300 MeV, the A resonance is often included as an explicit degree of freedom in the theory. 
The resulting chiral effective theory (xEFT) includes pion, nucleon, and A fields. When including the A as 
an explicit degree of freedom in the chiral Lagrangian, the counting scheme has to specify how the expansion 
parameter e = m^/K x sB is counted relative to 5 = A/A xSB . In the small scale expansion (SSE) HHem971L 
also called e-expansion, the pion mass and the Ma — mass difference are counted on the same footing, 
i.e. e ~ 8. The recently developed ^-expansion scheme, see [|Pas07l for a review and applications, counts the 
pion mass as e ~ 5 2 , which is the closest integer power relation between these parameters in the real world. 

Early calculations of the nucleon e.m. FFs in the SSE at order e 3 have been performed in Ref. [Ber98 ]. Be- 
cause such an approach is based on a heavy baryon expansion it is limited to Q 2 values much below 0.2 GeV 2 . 

Subsequently, several calculations of the nucleon e.m. FFs have been performed in manifestly Lorentz 
invariant x?T. Kubis and Meissner [KubOlJ performed a calculation to fourth order in relativistic baryon ;\PT, 
employing the infrared regularization (IR) scheme. They showed that the convergence of the chiral expansion 
is improved as compared to a heavy baryon ^PT results. Schindler et al. HSchi05ll also performed a manifestly 
Lorentz invariant calculation to fourth order, employing the extended on-mass-shell (EOMS) renormalization 
scheme. Both groups found that when including pion and nucleon degrees of freedom alone, one is not able 
to describe the nucleon e.m. FFs over a significant range of Q 2 . In both calculations, the proton electric FF 
would cross zero for Q 2 values as low as 0.4 GeV 2 . Both calculations also confirm that a realistic description 
of the nucleon e.m. FFs is only obtained once the vector mesons are included as explicit degrees of freedom 
in the chiral Lagrangian. The vector meson loop diagrams were found to play only a minor role, the dominant 
contribution coming from the pole diagrams, confirming the findings of VMD models and dispersion theory. 

The corresponding results for the nucleon e.m. FFs in both the IR and EOMS schemes are shown in Fig. [35] 
The covariant baryon ^PT results including vector mesons of Ref. UKubOll are shown in Fig. [35] at both third 
and fourth order. The electric FFs of proton and neutron require fixing one LEC for each, corresponding to 
the charge radii. One sees from Fig. [35] that the resulting fourth order results, including vector mesons, give 
a reasonably good description of the Q 2 dependence of the data up to Q 2 around 0.4 GeV 2 . For the magnetic 
FFs, at third order the two LECs are fixed from the corresponding proton and neutron magnetic moments, 
whereas at fourth order two more LECs are fixed from the magnetic radii. Also for the magnetic FFs, a good 
description is only obtained once the vector mesons are included. The results of Schindler et al. HSchi05l in a 
covariant %PT calculation and using a consistent power counting scheme which includes to fourth order both 
vector meson pole and loop contributions are also shown in Fig. [35] Again the explicit vector meson contri- 
butions play a major role at the higher end of this momentum transfer range in order to obtain a reasonable 
description of the data, as can be seen from Fig. [35] (compare dotted and dashed-dotted curves). 

4.5 Lattice QCD and chiral extrapolation 
4.5.1 Lattice simulations 

Lattice QCD calculations of nucleon structure quantities have matured considerably in the recent past. They 
provide an ab initio calculation of quantities such as the nucleon e.m. FFs from the underlying theory of QCD. 

Lattice QCD is a discretized version of QCD formulated in terms of path integrals on a space-time lat- 
tice HWil74B with only parameters the bare quark masses and the coupling constant. One recovers the con- 
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Figure 35: The nucleon e.m. FFs in the relativistic baryon xPT of [KubOl ] (IR scheme) and HSchi051 (EOMS scheme). 
The results of UKubOll including vector mesons are shown to third (dashed curves) and fourth (solid curves) orders. The 
results of [Schi05 ] to fourth order are displayed both without vector mesons (dotted curves) and when including vector 
mesons (dashed-dotted curves). For references to the data : Gm p (see Fig.©; Ge p (see Fig. [2 green triangles are data 
of [Cra06]); Gun (see Fig. ED); Geu (see Fig. |T8l and constraint from rms radius is given by green slopes [Kop97]). 



tinuum theory by extrapolating results obtained at finite lattice spacing a to a = 0. In order to perform the 
continuum extrapolation a separate calculation at several values of a is required. As lattice calculations neces- 
sarily are performed for a finite lattice size, one must keep the size of the box large enough to fit the hadrons 
inside the box. This requires to increase the number of sites as one decreases a. On the other hand, to keep 
finite volume effects small one must have a box that is much larger than the Compton wavelength of the pion. 
Present lattice QCD calculations take Lin^ ~ 5 where L is the spatial length of the box and the pion mass. 
As the computational costs of such calculations increase like m~ 9 , one uses quark mass values for the u and 
d quarks which are larger than in the real world. This enables the inversion of the fermionic matrix, which is 
needed for the calculation of hadronic matrix elements, with currently available resources. 

State-of-the-art lattice calculations for nucleon structure studies use a ~ 0.1 fm and L ~ 3 fm and reach 
pion mass values down to about 350 MeV. To connect those results with the physical world requires an extrap- 
olation down to the physical quark masses (note that m q is proportional to for small quark mass values). 
This so-called chiral extrapolation will be discussed further on. It is only very recently that pion mass values 
below 350 MeV [BerOO, Far04j| have been reached. This continuous effort is important to eliminate one source 
of systematic error associated with the extrapolation to the light quark masses. 

The bare coupling constant and quark masses are tuned as a changes to leave physical quantities un- 
changed. In a typical lattice calculation one starts by choosing the bare coupling constant g, which fixes the 
lattice spacing, and the bare masses for the u-, d- and s-quarks. One then computes a physical quantity such 
as the mass of the pion and the nucleon in lattice units as a function of the quark mass. The pion mass is used 
to fix the u- and d- quark masses (assumed degenerate) and the mass of the kaon or </> to fix the strange quark 
mass whereas the lattice spacing is determined by extrapolating the results, for instance, for the nucleon mass 
to the physical pion mass. Any other physical quantity in the light quark sector then follows. 

In the following, we will discuss lattice calculations for the (space-like) nucleon e.m. FFs. These calcula- 
tions require the evaluation of three-point functions, which involve two topologically different contributions as 
illustrated in Fig. [36] In the connected diagram contribution (left panel of Fig.l36l). the photon couples to one 
of the quarks connected to either the initial or final nucleon. The quark lines in Fig. [36] are understood to be 
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Figure 36: Diagrams illustrating the two topologically different contributions when calculating nucleon e.m. FFs in 
lattice QCD. Left (right) panels show the connected (disconnected) diagrams. Figure from [Boi06]. 

dressed with an arbitrary number of gluons exchanged between the quarks. If the fluctuations of such gluons 
into qq pairs are neglected, one speaks of the quenched approximation. The full QCD (unquenched) results 
include as well these sea-quark loop insertions into the gluon lines. The disconnected diagram (right panel 
of Fig. [36]) involves a coupling to a qq loop, which then interacts with the nucleon through gluon exchange. 
The disconnected diagram, which requires a numerically more intensive calculation, is at present neglected in 
most lattice studies. When taking the difference between proton and neutron e.m. FFs, i.e. for the isovector 
combination of nucleon e.m. FFs, the disconnected contribution drops out. Therefore, all following calcula- 
tions in which the disconnected diagram is neglected are applicable only to the isovector e.m. FFs. To directly 
calculate the proton and neutron e.m. FFs, involves the evaluation of the disconnected contribution, which 
awaits the next generation of dynamical-fermion lattice QCD simulations. 

The calculation of the connected diagram contribution to the nucleon e.m. FFs involves the computation of 
a sequential propagator, which can be done in two different ways. In an early pioneering work^of Leinweber, 
Woloshyn and Draper HLei91L this was done in the so-called fixed current approach, which requires the current 
to have a fixed direction and to carry a fixed momentum. This method allows one to use different initial or final 
states without requiring further inversions, which is the time-consuming part of the calculation. For a recent 
calculation of charge radii and magnetic moments of the whole baryon octet using this method, see [Boi06|. 
The drawback of this method is that a new calculation is required for each momentum transfer. 

More recently, a second method has been used by different groups to evaluate nucleon e.m. FFs, in which 
one fixes the initial and final states to have the quantum numbers of the nucleon. In this so-called fixed sink 
method, the current can couple to a quark line at any intermediate time slice, see Fig. [361 (left panel), carrying 
any possible value of the lattice momentum, which makes it the method of choice for a detailed study of the 
momentum transfer dependence of the nucleon e.m. FFs. 

The Nicosia-MIT group [Ale06a, Ale06bJ has performed a high-statistics calculation of nucleon isovector 
e.m. FFs in the fixed sink method, both in the quenched approximation and in full QCD, using two dynamical 
Wilson fermions, and for one value of the lattice spacing a, around 0.09 (0.08) fm for the quenched (un- 
quenched) results. The finite box size of length L imposes a smallest available non-zero momentum transfer, 
which for the quenched calculation is around Q 2 ~ 0.17 GeV 2 . The largest Q 2 value accessible is around 
Q 2 ~ 2 GeV 2 . Beyond such value, the Fourier transforms needed to evaluate the two- and three-point func- 
tions become noise dominated. Furthermore, the quenched calculation was performed for pion masses in the 
range m w ~ 410 — 560 MeV. The unquenched calculation was performed in the range of m n ~ 380 — 690 MeV. 

The lattice QCD results of [Ale06a] for the nucleon Dirac and Pauli isovector FFs are shown in Fig. [37] 

6 The first lattice QCD calculations for the pion e.m. FF were performed by Wilcox and Woloshyn [Wil85|, whereas the first 
attempt at a lattice QCD calculation for the proton electric FF was reported by Martinelli and Sachrajda |]Mar89l . 
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Figure 37: Lattice QCD results (from the Nicosia-MIT group [Ale06a]) for the isovector FFs (upper left panel) and 
F2 (lower left panel) as a function of Q 2 . Both the quenched results (Np = 0) and unquenched lattice results with two 
dynamical Wilson fermions (Np = 2) are shown for three different pion mass values. The right panels show the results 
for G\ (upper right panel) and Gj, (lower right panel), divided by the standard dipole FF, as a function of Q 2 in the 
chiral limit (using a linear extrapolation in m 2 ). The filled triangles show the experimental results for the isovector FFs 
extracted by interpolating the experimental data for the proton and neutron e.m. FFs. Figure from MAle06a11 . 



One observes that both the quenched and unquenched results for F± show only a very weak quark mass depen- 
dence in the range ~ 400 — 700 MeV. When comparing with experiment, one sees that both the quenched 
and unquenched lattice results of [Ale06a] largely overestimate the data for F^ . For F 2 y , one observes a 
stronger quark mass dependence, bringing the lattice results closer to experiment when decreasing m n . 

The two main uncertainties in this calculation are the continuum extrapolation (i.e. finite a effects) and 
whether one is close enough to the chiral limit (i.e. extrapolation in quark mass or m^). To check the latter, 
and to extrapolate the lattice results down to the physical pion mass value in order to directly compare with 
experiment, the Nicosia-MIT group uses a linear fit in (corresponding to a linear fit in the quark mass). 
Such a linear fit, which is supported by the lattice results in the range m n ~ 400 — 700 MeV, is also shown in 
Fig. [371 The thus extrapolated lattice results for F 2 y and G\ f are in agreement with experiment for Q 2 larger 
than about 0.3 GeV 2 . At smaller Q 2 , an agreement can also be expected as one can calculate in this range the 
m n dependence using %PT. We will discuss in the following that the pion loops lead to non-analytic behaviors 
in the quark mass, yielding e.g. a more rapid (than linear in m 2 ) variation of the isovector magnetic moment 
as one approaches the chiral limit. However, the linearly in m 2 extrapolated results for F^ still show strong 
disagreement with the data. This translates into an electric FF G V E which drops less fast than the dipole FF Gp>, 
whereas the data tell us that G V E drops faster than the dipole. It is puzzling that this strong disagreement is seen 
at larger values of Q 2 , where effects of pion loops are already suppressed, making it unlikely that the chiral 
extrapolation alone can explain this discrepancy. As both quenched and unquenched calculations of [Ale06aJ 
were only performed at one value of a, it would be very worthwhile, in order to shed light on this puzzle, to 
repeat such calculations for different values of a and check the continuum extrapolation. 
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Unquenched lattice calculations using two mass degenerate flavors of dynamical Wilson fermions have 
also been reported by the QCDSF Coll. [Goc06|. These results improve on previous calculations by the 
QCDSF Coll. HGoc05H which were performed using Wilson fermions in the quenched approximation. The Q 2 
dependence of the lattice results for the nucleon isovector FFs was parameterized (as a first approximation) 
in terms of a dipole behavior. For F± , the unquenched QCDSF lattice results find that the corresponding 
dipole mass becomes smaller with decreasing pion mass values. However at the smallest available pion mass 
values of around 340 MeV, the dipole mass reaches a value around 1.3 GeV, which still lies significantly above 
the experimental value of 0.843 GeV. This yields an isovector FF which has a too flat Q 2 dependence, 
confirming the puzzling finding of [Ale06aJ, which was also obtained for two dynamical Wilson fermions. 

Recently, unquenched lattice QCD calculations for the nucleon e.m. FFs have also been performed by 
the LHPC Coll. [Edw06] based on the Asqtad improved action, using different fermions for valence and sea 
quarks. This hybrid action uses for the valence quarks domain wall fermions, preserving chiral symmetry on 
the lattice. For the sea quarks, the configurations generated by the MILC Coll. [BerOOj are used, with two 
degenerate light and one strange staggered quarks, allowing for economical calculations. The discretization 
errors in this action are of order a 2 , in comparison with order a for the above discussed Wilson action. Although 
this action has generated quite a number of encouraging results when applied to nucleon structure studies, such 
as e.g. moments of unpolarized, helicity, and transversity distributions, see HEdw06ll for a recent overview, 
some controversy remains around the fourth root of the fermion determinant [Sha06 |. 

In Fig. [38l we show the unquenched lattice QCD results from the LHPC Coll. for the nucleon e.m. FFs, 
performed for one lattice spacing of a ~ 0. 125 fm, and for pion mass values in the range = 360 — 775 MeV. 
It is seen that in contrast to the above discussed Wilson results, this action yields a noticeable dependence on 

for the Dirac isovector FF F± at larger values of Q 2 . The Q 2 dependence of F^ at the smallest 
value of around 360 MeV is found to be in qualitative agreement with the data. One also sees from Fig. [38] 
that the isovector ratio F% /F± approaches the experimental result when decreasing m n . So far this is the 
only lattice calculation which yields a qualitative consistent picture for both F± and F 2 y . Evidently, it will 
be very worthwhile to corroborate the results at the lowest pion masses and improve their statistics in future 
calculations. If confirmed by higher statistics results, it remains to be understood why different actions may 
yield significantly different results, in particular for F^ . Unquenched calculations at a couple of different 
lattice spacings using different actions would be very helpful in this respect. 

4.5.2 Chiral extrapolations 

Present lattice calculations are possible for larger than physical quark masses, and therefore necessitate an 
extrapolation procedure in order to make contact with experiment. The extrapolation in the quark mass m q is 
not straightforward, because the non-analytic dependencies, such as ^ffrT q and h\m q , become important as one 
approaches the small physical value of m q . Therefore naive extrapolations often fail, while spectacular non- 
analytic effects are found in a number of different quantities, such as nucleon magnetic moments and charge 
radii, see e.g. Refs. ||Lei01[ |Hem02|. The %EFT, discussed in the Sect. 14.41 provides a framework to compute 
these non-analytic dependencies, for small quark masses. 

As an example, in the SSE (er-expansion) to order e 3 , the j*NN vertex has been calculated in HGoc05l 
through pion one-loop diagrams. Due to the pion loops, the isovector Dirac radius {r 2 ) \ acquires non-analytic 
dependencies in the quark mass. Its leading dependence in the pion mass is given by [Gas88, lLei93l : 



(r 2 )X = (r 2 ) lp -(r 2 ) ln = aM 



(4vr/.) 2 



f In ( + 0{ml) + ttA loops 




(44) 
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Figure 38: Lattice QCD results (from the LHPC Coll. Cdw0610 for the nucleon isovector FFs F± (left panel) and for 
the ratio F^f /F^ (right panel). The unquenched results using a hybrid action of domain wall valence quarks and 2+1 
flavor staggered sea quarks are shown for different values of and are compared with experiment (solid curve, using 
the parametrization of [Kel04]). Figure from IEd w06B . 



where the logarithmic term in m w is the leading non-analytic (LNA) dependence originating from the irN 
loop diagrams, and depends only on the nucleon axial coupling g& = 1.2695 and the pion decay constant 
f n = 92.4 MeV. Furthermore in Eq. (|44l) . fj, is the renormalization scale, and a ([i) is a LEC (evaluated at scale 
//). Note that, in contrast to most chiral extrapolations which contain finite terms of the form m 2 lnm^, the 
isovector radius diverges like logrn^, rendering the variation of the radius quite substantial near the physical 
pion mass. 

Analogously, the ttN loop diagrams give rise to non-analytic terms in the quark-mass expansion of the 
nucleon isovector magnetic moment. Its leading dependence in the pion mass is given by [Gas88J : 



K v = k p - K n = Kq - m 7v + 0( m i) + loops, (45) 

Y±K Jit 



where the LEC Kq corresponds to the isovector anomalous magnetic moment in the chiral limit. 

The LNA behavior of the isovector Pauli radius (r 2 )\ due to nN loops shows a l/m„ divergence in the 
chiral limit HGas88H : 

First attempts have been made to compare the lattice results for the pion mass dependence of the isovector 
magnetic moment, and the isovector Dirac and Pauli squared radii with the %PT results of HGoc06H in the SSE 
to order e 3 . For k v , allowing for one extra higher order parameter beyond e 3 , a four parameter ;vPT fit was 
performed in HGoc05l to the lattice data at relatively large values. It was found that for both k v and the 
isovector Pauli radius, (r 2 )^, the quenched and unquenched results from the Nicosia-MIT group [Ale0 6a J and 
the unquenched results from the QCDSF Coll. HGoc06l . at the smallest available values of indeed seem to 
follow the strong rise predicted by the %PT fit. On the other hand, the lattice results of [Ale06aJ and [Goc06J 
for the isovector Dirac radius, (r 2 )\ , show no strong indication of the logarithmic In divergence. The SSE 
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results of HGoc051l are not able to account for the lattice results for {r 2 )\ . 

One may of course wonder if any agreement or disagreement with %PT for m n values as large as 0.5 - 
1 GeV is very meaningful. Surely at such large m w values, higher order contributions not accounted for in e.g. 
the e 3 calculations are important^. A conservative strategy is to restrict xPT to its limited range of applicability 
and await lattice results for m w value below 300 MeV where the effect of higher order terms is still relatively 
small. Alternatively, one may choose to build upon xEFT and extend its range of applicability - leaving the 
domain of power counting - by resumming higher order effects using additional physics principles. 

One such strategy has been adopted in Refs. [Pas04, Hol05] by using analyticity to resum higher order 
(analytic) terms in m 2 to the nucleon magnetic moments. By requiring the anomalous magnetic moments to 
satisfy (a generalization of) the Gerasimov-Drell-Hearn sum rule [Ger66, Dre66J, a relativistic one-loop ttN 
calculation has the correct chiral behavior at the small values of m^, and yields a convergent 1/m 2 behavior 
at larger values of m n . It was found to yield a much smoother behavior for the magnetic moment than a 
truncated %PT calculation, while encompassing the correct behavior for small values. 

Another strategy has been pursued by the Adelaide group by modifying the one-loop ^PT results and taking 
into account the finite size of the nucleon through a finite range regularization procedure. This method was 
found successful when applying it to the calculation of the m w dependence of nucleon and A(1232) masses, 
see e.g. ULeiOOl ILei04ll . For (r 2 )\ , one may try in this spirit a modification of the %PT formula of Eq. (|44|) 



where A is a phenomenological cut-off which reflects the finite size of the nucleon. Such a fit for the isovector 
Dirac radius is shown in Fig.[39]and compared with the most recent unquenched lattice results using the hybrid 
action (domain wall valence quarks on top of a 2+1 flavor staggered sea) of the LHPC Coll. One firstly sees, 
that these lattice results do show appreciable m 2 variation over the pion mass range = 360 — 775 MeV and 
provide a first clear hint of the logarithmic m n divergence. As the pion mass approaches the physical value, 
the calculated nucleon size increases and approaches the correct value. Using the simple extrapolation formula 
of Eq. (147b . which has the In m w divergence at low m w values built in, one obtains a consistent description of 
the dependence of the lattice results using A ~ 500 MeV. 

Finally we would like to emphasize that presently there is no systematic framework for extrapolating lattice 
QCD results for FFs at values of Q 2 larger than about 0.3 GeV 2 , i.e. beyond the region where a ^PT expansion 
is expected to be applicable. The development of such a framework remains a challenge for future work. Even 
when lattice results become available for values below 300 MeV, at larger Q 2 , one is confronted with the 
problem of performing a chiral extrapolation (in the small scale m^) in the presence of a large scale Q 2 . A first 
attempt in this direction has been performed in [Mat05J, within the context of a light-front cloudy bag model. 

4.6 Generalized parton distributions (GPDs) 

So far we have discussed the N — > N transition as revealed with the help of the electromagnetic probe. 
By measuring the response of the hadron to a virtual photon, one measures the matrix element of a well- 
defined quark-gluon operator (in this case the vector operator q^q) over the hadronic state. This matrix 

7 See e.g. Ref. |McG06 |, where it was shown that the surprisingly good agreement of fourth-order ^PT when extrapolating lattice 
data for the nucleon mass out to large pion mass values (in the range 0.5 - 1 GeV) is spoiled once the fifth-order terms (due to 2-loop 
ttN diagrams) are included. 



as HDun02H : 




(47) 
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Dirac isovector charge radius 




Figure 39: Chiral extrapolation of the nucleon isovector Dirac radius (r 2 )\ . The unquenched results are from the LHPC 
Coll. IEd w06n . The experimental value is shown by the star. The curves are fits using the chiral extrapolation formula 
Eq. <|47]> of HDun02i Figure from HEdw06i 



element can be parametrized in terms of the nucleon e.m. FFs, revealing the quark-gluon structure of the 
nucleon. We are however not limited in nature to probes such as photons (or W, Z bosons for the axial 
transition). The phenomenon of asymptotic freedom of QCD, meaning that at short distances the interactions 
between quarks and gluons become weak, provides us with more sophisticated QCD operators to explore the 
structure of hadrons. Such operators can be accessed by selecting a small size configuration of quarks and 
gluons, provided by a hard reaction, such as deep inelastic scattering (DIS), or hard exclusive reactions such 
as deeply virtual Compton scattering (DVCS). We will be mostly interested here in DVCS reactions which 
are of the type 7*(<27i) + N(p) — > 'j(q') + N(p'), where the virtual photon momentum q^ is the hard scale. 
The common important feature of such hard reactions is the possibility to separate clearly the perturbative and 
nonperturbative stages of the interactions : this is the so-called factorization property. 

The all-order factorization theorem for the DVCS process on the nucleon has been proven in [Ji98al lCol99[ 
Rad98 |. Qualitatively one can say that the hard reactions allow one to perform a "microsurgery" of a nucleon 
by removing in a controlled way a quark of one flavor and spin and implanting instead another quark in the 
final nucleon. It is illustrated in Fig. [40] for the case of the DVCS process. The non-perturbative stage of 
such hard exclusive electroproduction processes is described by universal objects, so-called generalized parton 
distributions (GPDs) HMul941lJi97llRad96L see 0T98bl IGoeO 1 1 IDie03 1 lBel05l lJi04l for reviews and references. 

The nucleon structure information entering the nucleon DVCS process, can be parametrized at leading 
twist-2 level, in terms of four quark chirality conserving GPDs. The GPDs depend on three variables: the 
quark longitudinal momentum fractions x and £, and the momentum transfer Q 2 = —q 2 to the nucleon. The 
light-cone momentum fraction x is defined by k + = xP + , where k is the quark loop momentum and P 
is the average nucleon momentum P = (p + p ')/2, where p(p ') are the initial (final) nucleon four-momenta 
respectively, see Fig.|40l The skewedness variable £ is defined by q + = — 2£ P + , where q = p '— p is the overall 
momentum transfer in the process, and where 2£ — > x B /(l — x B /2) in the Bjorken limit: xb = Q\j (2p ■ qn) 
is the usual Bjorken scaling variable, with Q\ = —q\ > the virtuality of the hard photon. 

The DVCS process corresponds to the kinematics Q\^> Q 2 , M 2 , so that at twist-2 level, terms proportional 
to Q 2 /Ql or M 2 /Ql are neglected in the amplitude. In a frame where the virtual photon momentum q^ and the 
average nucleon momentum P^ are collinear along the z-axis and in opposite directions, one can parameterize 
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Figure 40: The "handbag" diagram for the nucleon DVCS process. Provided the virtuality of the initial photon (with 
momentum q^) is sufficiently large, the QCD factorization theorem allows to express the total amplitude as the convo- 
lution of a Compton process at the quark level and a non-perturbative amplitude parameterized in terms of generalized 
parton distributions (lower blob). The diagram with the photon lines crossed is also understood. 

the non-perturbative object entering the nucleon DVCS process as (following Ji HJi97ll fl 

i- J dy-e™ p+ y- (N(p')m-y/2) 1 -n^(y/2)\N(p))\ y+= , ±=0 

= H q (x,^Q 2 )N(p') 1 -nN(p) + E q (x,^Q 2 ) N(p) ta^^ N(p), (48) 

where ip is the quark field of flavor q, N the nucleon spinor, and n M is a light-cone vector along the negative 
^-direction. The Ihs of Eq. (|48l) can be interpreted as a Fourier integral along the light-cone distance y~ of 
a quark-quark correlation function, representing the process where a quark is taken out of the initial nucleon 
(having momentum p) at the space-time point y/2, and is put back in the final nucleon (having momentum p ') 
at the space-time point —y/2. This process takes place at equal light-cone time (y + = 0) and at zero transverse 
separation (y± = 0) between the quarks. The resulting one-dimensional Fourier integral along the light-cone 
distance y~ is with respect to the quark light-cone momentum xP + . The rhs of Eq. (|48T) parametrizes this 
non-perturbative object in terms of the GPDs H q and E q for a quark of flavor q. The quark vector operator 
(7 • n) corresponds at the nucleon side to a vector transition (parametrized by the function H q ) and a tensor 
transition (parametrized by the function E q ). Analogously, there are two GPDs corresponding to a quark axial 
vector operator (7 • 7175), which are commonly denoted by the polarized GPDs H q and E q . 

The variable x in the GPDs runs from —1 to 1. Therefore, the momentum fractions of the active quarks 
(x + £) for the initial quark and (x — for the final quark can either be positive or negative. Since positive 
(negative) momentum fractions correspond to quarks (antiquarks), it has been noted in [Rad96| that in this 
way, one can identify two regions for the GPDs: when x > £ both partons represent quarks, whereas for 
x < — £ both partons represent antiquarks. In these regions, the GPDs are the generalizations of the usual 
parton distributions from DIS. Actually, in the forward direction, the GPD H reduces to the quark (anti-quark) 
density distribution q(x) (q(x)) obtained from DIS: H q (x, 0, 0) = q(x), for x > 0; H q (x, 0, 0) = — q(— x), for 
x < 0. The GPD E is not measurable through DIS because the associated tensor in Eq. (|48l) vanishes in the 
forward limit (q — > 0). Therefore, E is a new leading twist function, which is accessible by measuring hard 
exclusive electroproduction reactions, such as DVCS. 

Besides coinciding with the quark distributions at vanishing momentum transfer, the GPDs have interesting 
links with other nucleon structure quantities. The first moments of the GPDs are related to the elastic FFs of 
the nucleon through model independent sum rules. By integrating Eq. (l48l over x, one obtains for any £ the 

8 In all non-local expressions we always assume the gauge link: Pexp(ig J dx^A^), ensuring the color gauge invariance. 
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following relations for a particular quark flavor [Ji97| 



f dxH q (x,£,Q 2 ) = F?(Q 2 ), J +1 dxE q (x,£,Q 2 ) = F|(Q 2 ) , (49) 

where Ff (F 2 5 ) represents the elastic Dirac (Pauli) FFs for the quark flavor q in the nucleon. These quark FFs 
are expressed, using SU (2) isospin, as flavor combinations of the proton and neutron elastic FFs as: 

F? = 2F lp + F ln + , F d = 2F ln + F lp + F( , (50) 

where Ff is the strangeness FF of the nucleon (which is neglected in the calculations discussed below). Rela- 
tions similar to Eq. (|5Q|) hold for the Pauli FFs F| . At Q 2 = 0, the normalizations of the Dirac FFs are given 
by: -^"(0) = 2 (Ff(Q) = 1) so as to yield the normalization of 2 (1) for the u (rf)-quark distributions in the 
proton. The normalizations of the Pauli FF at Q 2 = are given by F$(Q) = n q (for q = u, d), where k u , n d 
can be expressed in terms of the proton (k p ) and neutron (k ti ) anomalous magnetic moments as: 

k u = 1k v + n n = +1.673, K d = k p + 2k u = -2.033. (51) 

The above sum rules allow us to make a prediction for the nucleon e.m. FFs provided we have a model for 
the nucleon GPDs. Note that the sum rules of Eq. (|49l only involve valence quark GPDs, since the sea-quark 
and anti-quark contributions cancel each other in the sum rules. Since the results of the integration in Eq. (|49l 
do not depend on the skewness £0, one can choose £ = in these sum rules. We therefore only discuss the 
GPDs H and E at £ = in the following. 

In HDie05[ IGui051l . parameterizations of GPDs were developed which have a Regge behavior at small x 
and Q 2 , and which were modified to larger Q 2 behavior so as to lead to the observed power behavior of the 
FFs HBur03[|Bur04i A modified Regge parameterization for H and E was proposed in HGui05l : 

H«(x,0,Q 2 ) =q v {x)x a '^^\ E«(x,0,Q 2 ) = ^ (1 - x)^ q v {x) x a ' , (52) 

depending on 3 parameters. The Regge slope a' is determined from the Dirac radius, and two parameters rf 1 
and r] d , entering the GPD E, ensure that the x ~ 1 limit of E q has extra powers of 1 — x compared to that 
of H q . This results in a proton helicity flip FF F 2 which has a faster power fall-off at large Q 2 than F\, as 
observed experimentally. Furthermore, in Eq. (l52l) . the normalization factors N u and N d are given by : 

N u = [ dx (1 - x)^ u v (x) , N d = [ dx (1 - xf d d v (x) , (53) 

Jo Jo 

and guarantee the normalization condition for the GPD E q . 

Diehl et al. [Die05J chose a more general functional form for E q at the expense of more free parameters. 
In the following, we discuss the 'minimal' model with 3 parameters, and refer the interested reader to [Die05] 
for a study of more general functional forms. The 3 free parameters in the resulting modified Regge ansatz are 
to be determined from a fit to the FF data. 

In Fig. EH the proton and neutron Sachs electric and magnetic FFs are shown. One observes that the 
3-parameter modified Regge model gives a rather good overall description of the available FF data for both 



9 This is the simplest example of a so-called polynomiality condition when calculating moments of GPDs. 
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Figure 41: GPD calculation of Gup and Gmh relative to the dipole Gd (upper panels), ratio of Ge p /Gmp (lower 
left panel), and Geu (lower right panel), according to Ref. |Gui051 . The curves are a 3 parameter modified Regge 
parameterization : a' = 1.105 GeV -2 , rf = 1.713 and rf = 0.566. Data for Gmp are from [Jan66] (open squares), 
HLit70ll (open circles), HBer711 (blue solid stars), MBar731 (green open stars), HAnd94B (red solid circles), [Sil93] (red 
solid squares), according to the recent re-analysis of Ref. [Bra02]. Data for the ratio Ge p /Gmp are from HGay01[ (blue 
open triangles), ]Gay02[ (red solid squares), MPun05l (blue solid circles), and BCra06H (green solid triangles). The data 



for Gmu are from [XuOO] (red solid circles), MXu03ll (red solid squares), [Ank98] (open triangles), [Kub02] (green open 
stars), [Lun93 ] (open squares), MRoc921 (solid triangles), and HBro051 (blue solid stars). The data for Geu are from 
double polarization experiments at MAMI HHer991 IUst99l IBec991 lRoh99l IGla05H (red solid circles), NIKHEF HPas99H 
(green solid triangle), and JLab [ZhuOl , Mad03, War04 ] (blue solid squares). 



proton and neutron in the whole Q 2 range, using as value for the Regge trajectory a' = 1.105 GeV -2 , and the 
following values for the coefficients governing the x — > 1 behavior of the E-type GPDs: rf 1 = 1.713 and if = 
0.566. Note that a value i] q = 2 corresponds to a 1/Q 2 asymptotic behavior of the ratio F^/F^ at large Q 2 . 
The modified Regge GPD parameterization allows to accurately describe the decreasing ratio of Ge p /Gm p 
with increasing Q 2 , and also leads to a zero for Ge p at a momentum transfer of Q 2 ~ 8 GeV 2 , which will be 
within the range covered by an upcoming JLab experiment [Bra04J. 



4.7 Perturbative QCD (pQCD) 

The nucleon e.m. FFs provide a famous test for perturbative QCD. Brodsky and Farrar derived scaling rules 
for dominant helicity amplitudes which are expected to be valid at sufficiently high momentum transfers 
Q 2 HBro75ll . A photon of sufficient high virtuality will see a nucleon consisting of three massless quarks 
moving collinear with the nucleon. When measuring an elastic nucleon FF, the final state consists again of 
three massless collinear quarks. In order for this (unlikely process) to happen, the large momentum of the 
virtual photon has to be transferred among the three quarks through two hard gluon exchanges as illustrated 
in Fig.|42l This hard scattering mechanism is generated by valence quark configurations with small transverse 
size and finite light-cone momentum fractions of the total hadron momentum carried by each valence quark. 
The hard amplitude can be written in a factorized form HChe77a[ IChe77b[ |Efr79[ |Lep80[, as a product of a 
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Figure 42: Perturbative QCD picture for the nucleon e.m. FFs. The highly virtual photon resolves the leading three- 
quark Fock states of the nucleon, described by a distribution amplitude. The large momentum is transferred between the 
quarks through two successive gluon exchanges (only one of several possible lowest-order diagrams is shown). 

perturbatively calculable hard scattering amplitude and two distribution amplitudes (DAs) describing how the 
large longitudinal momentum of the initial and final nucleons is shared between their constituents. Because 
each gluon in such hard scattering process carries a virtuality proportional to Q 2 , this leads to the pQCD 
prediction that the helicity conserving nucleon Dirac FF F 1 should fall as 1/Q 4 (modulo In Q 2 factors) at suf- 
ficiently high Q 2 . Processes such as in Fig. |42l where the interactions among the quarks proceed via gluon 
or photon exchange, both of which are vector interactions, conserve the quark helicity in the limit when the 
quark masses or off- shell effects can be neglected. In contrast to the helicity conserving FF F l5 the nucleon 
Pauli FF F 2 involves a helicity flip between the initial and final nucleons. Hence it requires one helicity flip at 
the quark level, which is suppressed at large Q 2 . Therefore, for collinear quarks, i.e. moving in a light-cone 
WF state with orbital angular momentum projection l z = (along the direction of the fast moving hadron), 
the asymptotic prediction for F 2 leads to a 1/Q 6 fall-off at high Q 2 . 

We can test how well the above pQCD scaling predictions for the nucleon e.m. FFs are satisfied at currently 
available momentum transfers, see Fig. [43l One firstly sees that which has been measured up to about 30 
GeV 2 , displays an approximate 1/Q 4 scaling above 10 GeV 2 . For the proton ratio F 2p /Fi p , the data up to 5.6 
GeV 2 show no sign of a 1/Q 2 behavior as predicted by pQCD. Instead, the data show that the ratio F 2p /F lp 
falls less fast than 1/Q 2 with increasing Q 2 . Belitsky, Ji, and Yuan [Bel03| investigated the assumption of 
quarks moving collinearly with the proton, underlying the pQCD prediction. It has been shown in [Bel03J 
that by including components in the nucleon light-cone WFs with quark orbital angular momentum projection 
l z = 1, one obtains the behavior F 2 jF\ — > \n 2 (Q 2 / A 2 ) / Q 2 at large Q 2 , with A anon-perturbativemass scale0. 
Choosing A around 0.3 GeV, Ref. [Bel03] noticed that the data for F 2p / Fi p support such double-logarithmic 
enhancement, as can be seen from Fig. [43] (right panel). The arguments of [Bel03] still rely on pQCD and it 
remains to be seen by forthcoming data at higher Q 2 if this prediction already starts in the few GeV 2 region. 

Although at high enough Q 2 , the pQCD scaling predictions should set in, the available data for the nucleon 
e.m. FFs show that one is still far away from this regime. Nesterenko and Radyushkin [Nes83J argued that 
the above described hard scattering mechanism is suppressed at accessible momentum transfers relative to 
the Feynman mechanism [ |Fey72[ , also called soft mechanism. The soft mechanism involves only one active 
quark, and the FF is obtained as an overlap of initial and final hadron WFs. The hard scattering mechanism on 
the other hand, involving three active quarks, requires the exchange of two gluons, each of which brings in a 
suppression factor a s /n ~ 0.1. One therefore expects the hard scattering mechanism for F\ p to be numerically 
suppressed by a factor 1/100 compared to the soft term, see also [Bo l95[|Kro91|| . Even though the soft mech- 
anism is suppressed asymptotically by a power of 1 /Q 2 relative to the hard scattering mechanism, it may well 



10 In llRal04l lBro03L it has also been discussed that inclusion of quark orbital angular momentum yields a ratio F2 P /Fi p which 
drops less fast than 1/Q 2 with increasing Q 2 . 
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Test of the scaling behavior of the proton FFs. Upper left panel : proton Dirac FF multiplied by Q 4 . Lower 
ratio of Pauli to Dirac proton FFs multiplied by Q 2 . Right panel : test of the modified scaling prediction for 
F2p/F\ p of BBeKBL The data for F\ p are from [Sil93] (solid squares). Data for the ratio F2 P / F\ p on both panels are 
from [Pun05] (blue solid circles), [GayOl] (red open triangles), and |Gay02] (red solid squares). The curves on the left 
panels represent the calculation based on the three parameter modified Regge GPD parametrization of [Gui05 ]. 



Figure 43 
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dominate at accessible values of Q 2 . In |Nes 83J. the soft contribution to the nucleon e.m. FFs was estimated 
using a model based on local quark-hadron duality, and was found to yield an approximate 1/Q 4 behavior in 
the range Q 2 ~ 10 — 20 GeV 2 , in qualitative agreement with the F lp data. 

In a more recent work, the soft contribution was evaluated within the light-cone sum rule (LCSR) ap- 
proach [Bra06] . Using asymptotic DAs for the nucleon, the LCSR approach yields values of G m p and Gmu 
which are within 20 % compatible with the data in the range Q 2 ~ 1 — 10 GeV 2 . The electric FFs however were 
found to be much more difficult to describe, with Geh overestimated, and Ge p /Gm p near constant when using 
an asymptotic nucleon DA. Only when including twist-3 and twist-4 nucleon DAs within a simple model, is a 
qualitative description of the electric proton and neutron FFs obtained. Such higher twist components hint at 
the importance of quark angular momentum components in the nucleon WF. 

In Sect. I4.6L we have shown that the nucleon e.m. FFs can be obtained from model independent GPD sum 
rules. These GPDs, represented by the lower blob in Fig.|40l are non-perturbative objects which include higher 
Fock components in the nucleon WFs. One can use a GPD parametrization to provide an estimate of the soft 
contributions, and expects this non-perturbative approach to be relevant in the low and intermediate Q 2 region 
for the FFs. This is shown in Fig. [43] (solid curves) from which one sees that the GPD Regge parametrization 
discussed above is able to explain at the same time an approximate 1/Q A behavior for F lp and a behavior for 
F 2p /F lp which falls less steep than l/Q 2 . Forthcoming experiments at the Jefferson Lab 12 GeV facility will 
extend the data for F 2p / F lp to Q 2 values around 13 GeV 2 . Such measurements will allow to quantify in detail 
the higher Fock components in the nucleon WF (which are all included in the nucleon GPD) versus the simple 
three-quark Fock component, and to map out the transition to the perturbative QCD regime. 
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5 Conclusions and outlook 



The increasingly common use of the double-polarization technique to measure the nucleon FFs, in the last 15 
years, has resulted in a dramatic improvement of the quality of all four nucleon e.m. FFs, Ge p , Gm p , Geu an d 
Gmti- It has a l so completely changed our understanding of the proton structure, having resulted in a distinctly 
different Q 2 - dependence for Ge p and Gm p , contradicting the prevailing wisdom of the 1990's based on cross 
section measurements and the Rosenbluth separation method, namely that Ge p and Gm p obey a "scaling" 
relation hGe p ~ Gm p - A direct consequence of the faster decrease of Ge p revealed by the JLab polarization 
experiments was the disappearance of the early scaling F 2 /Fi ~ l/Q 2 predicted by perturbative QCD. 

The main origin of this abrupt change in results is now understood in simple terms. The faster decrease 
of Ge p reduces its contribution to the cross section significantly below the natural ratio prevailing at small 
Q 2 , namely G 2 Ep /G 2 M ~ l//i 2 . At the highest Q 2 for which we now have polarization data, 5.6 GeV 2 , the 
contribution from the electric FF to the cross section is less than 1%. It has been realized in recent years, that 
to extract Ge p from Rosenbluth separations at larger Q 2 requires a much better quantitative understanding of 
several of the radiative corrections contributions, including in particular the one due two hard photon exchange. 
As discussed in section 1331 there are currently differences of order several %s between the results of various 
radiative correction calculations. The two hard photon correction by itself might explain the whole discrepancy 
between Rosenbluth and recoil polarization results, but it does not affect recoil polarization results measurably, 
because these are measurements of ratios of FFs and both FFs are, in first order, modified similarly. Until the 
origin of the difference between cross section- and polarization results is understood in full quantitative detail, 
it is safe to take the polarization results as the closest to the real, Born approximation, proton FFs. 

The use of the polarization technique has also resulted in a constant progress in the measurement of Geu, 
which is intrinsically more difficult to obtain because of the smallness of this form factor, due to the overall 
zero charge of the neutron. Recent times have seen the maximum Q 2 for which we have polarization FFs grow 
to 1.5 GeV 2 , with new data obtained and under analysis up to 3.4 GeV 2 , and several experiments planned or 
proposed to significantly higher Q 2 values. Important progress has been made for G Mn too, with new data 
with much improved error bars up to 4.8 GeV 2 . 

A basic understanding of the nucleon e.m. FFs can be obtained within the VMD picture, in which the pho- 
ton couples to the nucleon through the exchange of vector mesons. Dispersion analyses build on this picture 
by including, besides vector mesons, also non-resonant contributions in the coupling of the photon with the 
nucleon. The state-of-the-art dispersion analyses, which include the 2n continuum in the isovector spectral 
function, and the KK and pn continua as input in the isoscalar spectral function, are found to yield a reason- 
ably good overall description of the data for all four nucleon e.m. FFs using a 15 parameter fit. 

The effect of pionic degrees of freedom in the nucleon e.m. FFs can be systematically calculated within 
chiral effective field theory. The latest relativistic %EFT calculations found that calculations based on pions 
and nucleons alone are not able to explain the nucleon e.m. FFs. Only upon inclusion of explicit vector meson 
degrees of freedom, these calculations were found to describe the FFs in the range Q 2 < 0.4 GeV 2 . 

Calculations of the nucleon e.m. FFs within constituent quark models have highlighted the role of relativity 
when trying to arrive at a microscopic description of nucleon FFs based on quark degrees of freedom in the 
few GeV 2 region. Although a complete calculation is independent of the specific choice of relativistic form 
chosen to describe the dynamics, present approximations destroy this independence. 

We have also reviewed the recent progress made by lattice QCD calculations of the nucleon e.m. FFs. In 
present lattice simulations, disconnected diagrams, which are numerically more intensive, have not yet been 
evaluated. Therefore the current lattice calculations are for the isovector combination of nucleon FFs where 
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this disconnected diagram contribution drops out. They are performed at pion mass values above about 350 
MeV and have to be extrapolated to the physical pion mass to allow for a comparison with experiment. The 
calculations using dynamical Wilson quarks are able to provide a reasonably good description of the isovector 
Pauli FF F 2 V . At the lower Q 2 , the non-analytic terms in the chiral extrapolation were found to be important 
to arrive at a description of the isovector magnetic moment and Pauli radius. The present dynamical Wilson 
calculations largely overestimate the isovector Dirac FF F± however. In contrast, the results using a hybrid 
action, consisting of domain wall valence quarks and staggered sea quarks, are compatible for F 2 V , but differ 
from the Wilson results for F± . The F± results using the hybrid action are in qualitative better agreement with 
the data, and provide a first clear hint of the logarithmic m n divergence in the isovector Dirac radius. 

The quark structure of the N — > N electromagnetic transition can be accessed in hard scattering pro- 
cesses such as deeply virtual Compton scattering. The non-perturbative information in this process can be 
parameterized in terms of GPDs, and the nucleon e.m. FFs can be obtained as first moment sum rules in the 
quark longitudinal momentum fraction of these GPDs. In particular, the GPDs contain the information on 
the quark distribution of the nucleon anomalous magnetic moment which can not be accessed from inclusive 
deep inelastic scattering experiments. We discussed a GPD parameterization, and found that the first moment 
of such a 3-parameter modified Regge parameterization yields a good overall description of the nucleon e.m. 
FF data over the whole available Q 2 range. The modified Regge GPD parameterization predicts that, at mod- 
erately large Q 2 values, F lp follows an approximate l/Q 4 scaling, whereas F 2p /Fi p drops less fast than the 
1/Q 2 pQCD behavior, in agreement with the polarization data. It furthermore predicts that Ge p reaches a zero 
around Q 2 ~ 8 GeV 2 , which is in the reach of an upcoming experiment. 

We like to end this review by spelling out a few open issues and challenges in this field : 

1. Quantifying the two-photon exchange processes both experimentally and theoretically 

In order to use electron scattering as a precision tool, it is clearly worthwhile to arrive at a quantita- 
tive understanding of two-photon exchange processes. This calls for detailed experimental studies, and 
several new experiments are already planned. Differences between elastic e~ and e + scattering directly 
access the real part of the two-photon exchange amplitude. The predicted small effect of two-photon pro- 
cesses on the polarization data can be checked by measuring the e dependence in polarization transfer 
experiments. These upcoming experiments also call for further refinements on the theoretical side. 

2. Dispersion analyses 

In the present dispersion analysis for the nucleon e.m. FFs, the pQCD limit F 2 jF\ ~ l/Q 2 was built in 
as a constraint, although the data do not support this limit. It might therefore be worthwhile to investigate 
how the dispersion fits change by removing this bias from the analysis at the largest available Q 2 values. 

3. Relativistic constituent quark model calculations 

The quality of the commonly introduced impulse approximation when describing nucleon FFs in rela- 
tivistic CQMs may differ between different forms. As a next step for CQMs, it would clearly be worth- 
while to investigate the approximations made in the current operator within each form, and quantify e.g. 
the effect of explicit two-body currents. 

4. Lattice calculations and chiral extrapolations 

(a) One would clearly like to understand the present disagreement between the unquenched lattice 
predictions when using different actions (Wilson action vs. hybrid action). Understanding the 
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structure of the FFs at low Q 2 , such as the Dirac charge radius, depends crucially on the effect 
of pion loops, which yield strong non-analytic dependence (in particular a log m n singularity). 
At the larger Q 2 values, between 0.5 and 2 GeV 2 , the effects of pion loops are expected to be 
reduced however. One may therefore suspect that the differences between different actions at 
larger Q 2 values are due to discretization errors, which have only partly been studied. Unquenched 
calculations at a couple of different lattice spacings using both the Wilson and hybrid actions would 
be very helpful to shed a further light on this issue. 

(b) In order to provide predictions for proton and neutron e.m. FFs separately, the calculation of the 
disconnected diagrams awaits the next generation of dynamical fermion lattice QCD simulations. 

(c) A further challenge for the lattice calculations is a fully consistent treatment of both valence and 
sea quarks which respect chiral symmetry on the lattice. 

(d) As future lattice calculations for pion mass values around and below 300 MeV become available 
for FFs in the range Q 2 ~ 0.5 GeV 2 , i.e. beyond the range where present %PT calculations are 
applicable, one is confronted with a two-scale problem. A challenge is to theoretically study the 
extrapolation (in the small scale m n ) in the presence of a (moderately) large scale Q 2 . 

5. Precision measurements in the low Q 2 regime 

Precision measurements of the nucleon e.m. FFs in the Q 2 range below 0.5 GeV 2 , may bring the effects 
of the pion could sharper into focus. In this respect, new measurements for Geu are needed to better 
quantify the conjectured "bump" structure in G En around Q 2 ~ 0.3 GeV 2 . 

6. Extending the FF measurements to larger Q 2 

The anticipated upgrade of JLab to 12 GeV beam energy, offers promises of measurements of all four 
FFs up to or larger than ~10 GeV 2 . It is unlikely that we will see indications of a clear departure from 
the soft physics dominance to the hard collision regime of pQCD. However these data will constrain 
the parameterizations of GPDs and yield information on the spatial distribution of partons which carry a 
large momentum fraction of the nucleon momentum, i.e. partons with x ~ 1. 

The recent unexpected results in the nucleon e.m. FFs using double-polarization high-precision experi- 
ments, have challenged our theoretical understanding of the structure of the nucleon. They have triggered 
several new theoretical developments, which were reviewed in this work. As a result of the unexpected find- 
ings, several further experiments are planned, which will bring the quark-gluon structure of the most common 
constituent of visible matter in the universe into sharper focus. 
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